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In  an  effort  to  identify  the  dominant  source  of  the  radiated  aero-acoustic  noise  pro¬ 
duced  by  high-speed,  heated  jets  and  develop  control  strategies  to  reduce  it,  two  primary 
experiments  are  led  in  parallel  using  a  2  inch  diameter  nozzle  operating  at  Mach  0.6:  I) 
An  investigation  of  the  aeroacoustic  sources  of  two  ’acoustically  matched’  Mach  0.6  jets, 
at  temperature  ratio  of  7>=0.93  and  7Y= 1.7,  is  conducted  to  identify  distinctions  in  con¬ 
tribution  to  the  acoustic  far-field  spectrum.  The  fluctuating  pressure  sampled  near  the 
exit  of  the  jet(x/D=2,  r/D  ~1),  is  shown  to  sense  the  dominant  flow  features  intrinsic  to 
each  velocity  field.  POD  analysis  of  the  velocity  fields  reveals  the  presence  of  larger  scaled, 
more  energetic  structures  in  the  hot  jet,  as  well  as  a  dominant  Fourier-column  mode-like 
structure,  and  a  Fourier-helical  mode-like  structure  in  the  cold  jet.  Fourth  order  correla¬ 
tions  of  the  self  noise  contributions  are  shown  to  be  dominant  in  both  jets,  with  the  hot  jet 
exhibiting  a  greater  shear  noise  efficiency,  and  cold,  a  greater  self  noise  contribution.  The 
low  frequency  emission  of  the  shear  noise  terms  is  linked  to  the  increase  in  the  acoustic 
spectrum  of  the  hot  jet  at  these  frequency,  where  as  the  decrease  at  high  frequency  is 
characteristic  of  the  self  noise  weighting.  The  pronounced  effect  of  increased  low  frequency 
noise  is  not  deemed  likely  as  solely  resulting  from  the  changes  in  the  contribution  of  the 
shear  noise  sources,  and  as  such,  points  to  the  contribution  of  the  second  entropy  source 
term.  II)A  dual-time  PIV  investigation  of  the  sound  producing  region  (x/D= 3  thru  10), 
(7>=0.93),  is  performed  to  capture  the  instantaneous  Eulerian  acceleration  field.  An  iden¬ 
tification  method  is  implemented  to  solve  for  the  coefficients  of  the  ordinary  differential 
equation  (ODE)  derived  from  a  Galerkin  projection  of  the  Navier-Stokes  equations  onto  a 
basis  of  POD  eigenfunctions,  describing  the  evolution  of  the  flow.  An  extensive  database 
of  the  velocity  and  acceleration  fields  is  used  as  training  data  for  the  development  of  an 
empirical  low-order  dynamical  system  (LODS),  for  eventual  use  in  closed-loop  flow  control 
applications.  These  predictive  models  show  a  good  agreement  with  the  original  data  when 
the  truncation  of  the  system  is  minimized.  It  is  found  beneficial  to  retain  more  modes  in 
a  quadratic  model  than  to  increase  the  order  of  the  system,  which  restricts  the  number  of 
modes  that  can  be  used  in  the  modeling.  The  higher  modes  are  found  to  be  responsible  for 
much  of  the  higher  frequency  content  in  the  models,  which  make  them  more  representative 
of  the  actual  flow. 


I.  Introduction 

Since  its  advent,  and  throughout  it’s  life  cycle  from  gas  turbine  to  combustion,  the  jet  engine  has  been 
the  focus  of  numerous  research  efforts.  Of  particular  interest  is  the  mechanism  by  which  the  exhausted 
plume  generates  noise.  At  the  heart  of  this  problem  is  the  challenge  of  deciphering  what  events  in  these 
high-speed  flows  arc  most  efficient  at  radiating  highly  intense  acoustic  pressure  fluctuations,  that  translate 
to  the  far-field  as  broadband  noise.  The  noise  due  to  current  military  and  commercial  high-speed  jet  aircraft 
engines  can  reach  upwards  of  120  ( dB )  on  take  off  (depending  on  proximity),  with  the  exhausted  plume 
being  the  most  significant  contributor.2  With  regulations  pertaining  to  noise  from  commercial  jet  aircraft 
becoming  more  stringent,  efforts  to  mitigate  noise  levels  are  a  major  concern.  Military  applications  remain 
focused  on  improving  stealth  technology  on  a  range  of  aircraft,  from  unmanned  air  vehicles  (UAV)  to  fighter 
aircraft.  Where  radar  signature  is  now  minimized  on  the  most  recent  aircraft,  their  noise  signature  has  yet 
to  be  considerably  altered. 

The  noise  ’source’  in  jets  has  been  characterized  by  the  generation,  subsequent  interaction,  and  evolution 
of  the  coherent  turbulence  structure  within  the  shear  layer  produced  by  the  mixing  between  the  high- 
momentum,  high-temperature  core  jet  flow  and  the  lower-speed,  cold  entrained  flow,  as  cited  by  Ffowcs 
Williams  &  Kcmpton.3  Research  at  the  most  fundamental  level  involves  the  use  of  circular  nozzles  (coaxial, 
over,  ideally,  and  under  expanded)  and  primarily  cold  (ambient  temperature)  core  flow.  The  seemingly 
subtle  differences  between  cold  and  hot  jets  can  prove  to  be  extremely  complex  in  the  experimental  arena, 
as  solutions  to  Lighthill’s  wave  equation  for  noise  caused  by  unsteady  fluid  motion,  involves  making  certain 
assumptions  that  do  not  apply  at  the  elevated  temperatures  at  which  realistic  engine  cycles  occur.  Terms 
that  can  simply  be  deemed  negligible  under  the  assumption  of  being  small  in  the  cold  jet,  must  now  be 
included  through  measurement.  Unfortunately  for  the  researcher,  the  development  of  measurement  tools 
lags  behind  the  necessity  for  required  information.  We  must  therefore  rely  on  numerical  tools  such  as  Direct 
Numerical  Simulation  (DNS)  and  Large  Eddy  Simulation  (LES)  to  provide  answers.  These  tools  provide 
solutions  to  the  equations  of  motion,  but  at  times  leaves  much  to  be  desired  in  terms  of  the  physical  flow 
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dynamics,  as  full  3-D  simulations  are  costly  both  in  capital  and  time.  Many  studies  have  recently  moved 
towards  using  this  information  to  validate  experimental  results.  However,  in  the  case  of  gleaning  insight  into 
the  turbulent  mixing  mechanism  responsible  for  noise  generation  in  the  hot  jet,  the  coupling  of  numerical 
and  experimental  data  appears  to  be  the  only  current  solution. 

The  area  of  flow  control,  geared  toward  the  reduction  of  jet  noise,  is  in  its  early  stages  of  development. 
Currently  implemented  strategics  consist  of  passive  control  devices,  commonly  known  as  chevrons  or  tabs. 
These  serve  to  enhance  mixing  between  the  high-momentum,  high-temperature  core  jet  flow  and  the  lower- 
speed,  cold  entrained  flow.  These  devices  have  shown  some  promise  in  reducing  the  noise  radiated  by  these 
high-speed  jet  exhausts  without  compromising  the  thrust  efficiency,  notably,  reducing  low  frequency  noise  at 
aft  angles  while  increasing  high  frequency  noise  at  broadside  angles  relative  to  the  jet.  However,  significant 
noise  mitigation  has  yet  to  be  achieved. 

Flow  control  at  the  first  stages  of  the  design  process  is  needed  for  such  complex  problems  where  so  many 
disciplines  come  into  play;  from  propulsion  to  thermodynamics,  turbulence,  acoustics  and  control  theory. 
The  more  realizable  alternative  of  ’active’  flow  control  has  now  gained  considerable  focus.  To  date,  active 
flow  control  techniques  for  jet  noise  reduction  have  not  yet  provided  enough  consistency  and  cost-efficiency 
to  be  found  in  industrial  applications.  Development  of  theses  strategies  represents  the  solution  with  the 
most  potential  for  furthering  and  optimizing  the  noise  reduction  capability  provided  by  passive  techniques. 
Ultimately,  both  active  and  passive  flow  control  techniques  could  be  complementary  in  accomplishing  this 
challenging  goal. 


II.  Jet  Noise 

The  turbulent  motion  of  shear  flows  has  long  proved  challenging  to  researchers.  On  the  heels  of  the 
work  of  Sir  James  Lighthill  in  identifying  the  mechanism  by  which  sound  is  generated  acrodynamically,6 
and  Ffowc  Williams  &  Kcmpton3  reporting  that  jet  noise  production  was  biased  towards  the  contribution 
of  large-scale  turbulence,  many  groups  set  out  to  identify  orderly  patterns  in  the  flow  as  ’sources’  of  sound, 
in  hopes  of  controlling  it. 

The  establishment  of  Lighthill’s  aerodynamic  noise  theory6,7  stands  as  one  of  the  more  significant  contri¬ 
butions  to  identifying  the  role  of  shear  flows  in  generating  vortical  structures,  and  the  noise  caused  by  this 
unsteady  fluid  motion.  An  acoustic  approximation  to  momentum  transport,  Lighthill’s  theory  describes  the 
production  of  ’aero-acoustic’  noise  by  high-speed  jets  through  the  time  dependent  kinematics  which  generate 
pressure  fluctuations  as  they  evolve  downstream  of  the  jet  exit,  ultimately  radiating  to  the  far  field  as  noise. 
Solutions  to  Lighthill’s  wave  equation  realizes  the  far-ficld  fluctuating  quantities  (t.e.  acoustic  pressure, 
density),  but  involves  making  certain  assumptions  that  do  not  apply  at  the  elevated  temperatures  at  which 
realistic  engine  cycles  occur.  Where  the  axisymmetric  jet  is  of  concern,  it  is  the  factors  surrounding  the 
propagation  of  this  sound  to  the  far-ficld,  as  well  as  the  entropy  contribution  in  heated  jets,  that  remains 
not  clearly  identifiable. 


A.  The  acoustic  analogy 

Lighthill’s  aerodynamic  noise  theory6,7  identifies  a  relationship  between  the  sound  source  in  free  shear  flows 
and  the  resulting  distribution  of  external  pressurc/dcnsity  fluctuations.  An  analogy  is  drawn  between  the 
non-linear  equations  of  motion  and  the  linear  theory  of  acoustics  to  form  the  inhomogeneous  wave  equation, 


d2p 


~  a20V2p  = 


d2Txl 
doc  i  doc  j 


(1) 


where  the  left  hand  side  of  the  equation  represents  acoustic  wave  propagation,  and  on  the  right  the  ’source’ 
term,  Lighthill  stress  tensor  TtJ  is  defined  as, 

Tij  =  puxUj  +  {p  —  o,Qp)dij  —  T{j  (2) 
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In  the  simple  ease  of  a  cold,  subsonic  jet  the  stress  tensor  can  be  reduced  to  the  contribution  of  the  reynolds 
stresses  only, 


Txj  =  putUj 


(3) 


by  a  series  of  assumptions.  Firstly,  rxj  the  viscous  stress  tensor,  can  generally  be  neglected  as  jet  noise 
is  considered  to  be  driven  by  an  inviscid  process.24  Secondly,  for  mildly  heated  isothermal  jets  the  term 
commonly  referred  to  as  the  ’entropy’  term  (p-a^p)^  is  also  considered  small,  as  (p'  =  a0p')  is  an  accurate 
iscntropic  assumption  in  isothermal  jets,  and  neglected.74  The  solution  to  Lighthill’s  wave  equation  for 
far-ficld  pressure  fluctuation  can  be  expressed  in  free  space  as  the  volume  integral, 


p(x.i) 


JL  T)_*_ 

47 Tjvdxidxjyy'  ;|x-y| 


(4) 


for  the  contribution  resulting  from  V,  an  average  eddy  volume.  The  variables  x  and  y  are  location  of 
observation  reference  point,  and  fluid  volume  respectively,  and  t  —  t  —  is  the  retarded  time. 


Lillcy 19  in  an  attempt  to  include  flow  variables  associated  with  the  propagation  of  sound  through  the  flow 
field  (and  subsequent  interaction  of  sound  waves  and  flow,  flow-acoustic  interaction )  established  that  contrary 
to  Lighthill’s  original  formulation  involving  only  terms  of  quadrupole  nature,  there  existed  a  component  of 
the  pressure-density  relationship  more  appropriately  represented  by  dipole  sources.  The  term  (p  —  alp)6lJ 
argued  by  Lighthill  to  be  associated  with  non-isentropic  density  fluctuations  produced  by  viscous  effects  was 
shown  to  also  have  an  iscntropic  component,  resulting  in  an  important  dipole  noise  source  in  the  presence  of 
density/temperature  fluctuations.  Howe20  in  a  reformulation  with  emphasis  placed  on  this  ’excess’  noise  of 
hot  jets  found  terms  of  similar  dipole-like  contribution,  discussed  in  further  detail  later  in  the  next  section. 


B.  Hot  jets 

Most  effort  toward  the  identification  of  turbulent  sources  of  sound  and  development  of  jet  noise  theory 
has  been  focused  on  the  subsonic,  cold  jet  due  to  it  measurable  characteristics.  However,  at  the  elevated 
temperatures  at  which  realistic  engine  cycles  occur,  it  has  been  shown  that  noise  sources  are  altered,  whereby 
entropy  fluctuations  are  now  prominent,  producing  a  noticeably  different  far-ficld  spectrum. 

The  assumption  previously  made  in  Lighthill’s  acoustic  analogy  regarding  the  (p-a^p)6ij  term,  no  longer 
holds  with  the  addition  of  heat.  Lighthill  briefly  discussed  that  in  the  heated  jet,  this  term  could  be  replaced 
by  [1  —  (oq /a2)]  V2p,  and  quadrupole  type  sources  could  again  be  assumed,  however  this  formulation  was  later 
shown  to  be  incomplete.  In  attempts  to  adapt  LightilPs  original  theory  to  flows  with  mean  flow-acoustic 
interaction  (i.e.  sources  in  a  convecting  medium)  researchers  (Lilley, 19  Goldstein,17  Howe20)  derived  an 
inhomogeneous,  convective  wave  equation  from  the  Navier-Stokes  equations.  Mainly  varying  in  assumptions 
made  and  choice  of  independent  variables  in  the  final  approximation  to  far-field  pressure,  these  theories 
result  in  greater  complexity  involving  solutions  to  the  full  unsteady  equations  of  motions. 

Howe,20  in  a  reformulation  with  emphasis  placed  on  the  ’excess’  noise  arising  in  heated  flow,  used  the 
velocity  potential  and  stagnation  enthalpy  as  independent  variables  to  derive  an  analogy  relating  the  stagna¬ 
tion  enthalpy  to  vorticity  and  specific  entropy.  To  demonstrate  its  application,  a  model  problem  is  considered 
in  which  the  rate  of  heat  energy  absorption  by  a  fluid  particle  can  be  expressed  in  terms  of  the  temperature 
gradient  across  the  fluid, 


pT—  =div(/Cgradr)  (5) 

where  K  denotes  the  thermal  conductivity  of  the  fluid  and  S  specific  entropy.  Combining  the  continuity 
equations  and  (Eqn.  5),  and  through  a  series  of  assumptions  limiting  the  thermodynamic  quantities  to  first 
order,  the  acoustic  equation  can  now  be  written  as, 


( 


\_(P_ 
al  dt2 


d_ 

dt 


(6) 
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where  S'  represents  the  fluctuating  component  of  specific  entropy,  and  p  acoustic  pressure  perturbations. 
A  solution  for  the  radiated  pressure  field  resulting  from  this  expansion  and  contraction  of  the  fluid  is  then 
realized,  through  the  use  of  an  appropriate  Green’s  function. 

Similar  results  for  the  far-field  pressure  were  obtained  in  the  alternative  matched  expansions  analysis  of 
Landau  &;  Lifshitz.37  Although  a  simple  model,  the  contribution  of  entropy  fluctuations  to  far-field  sound 
is  highlighted.  Howe20  concluded  that  these  fluctuations  only  radiate  noise  in  the  presence  of  a  mean  flow 
pressure  gradient. 

Lush  et  al .36  also  showed  the  plausible  existence  of  an  additional  entropy  fluctuation  source  in  heated 
flows  by  application  of  the  isentropic  relationship,  S  =  cvln(p/p'|r)  -I-  const ,  for  a  perfect  gas,  where  S  is 
entropy  and  7  =  cp/cv.  From  this, 


dp 

dt 


1  dp 
a2  dt 


p_dS 
cv  dt 


The  time  derivative  of  the  (p  -  alp)  term  from  Light  hill’s  stress  tensor  then  becomes, 

\  /vi  ?  . 


(8) 

(9) 


where  a  =  (7 p/p)^  is  local  sound  speed.  The  derivation  highlights  the  importance  of  the  second  term  on  the 
right,  as  large  temperature  gradients  between  core  flow  and  ambient  entrained  flow  would  invariably  create 
large  entropy  fluctuations.  It  also  validates  the  assumption  that  allows  the  (p  —  a^p)  term  to  be  neglected  as 
small  in  cold  jets,  as  a  =  ao  and  entropy  is  generally  considered  uniform  (dp  —  aldp)  across  the  flow  field. 


C.  Far-field  noise 

Even  with  this  additional  ’source’  present,  the  overall  contribution  of  heat  to  the  far-field  noise  is  ultimately 
a  function  of  exit  velocity.  To  illustrate  this  effect,  the  acoustic  Mach  number  is  introduced,  Ma  =  U/a 0, 
where  U  is  the  jet  exit  velocity  and  ao  the  ambient  sound  speed.  At  fixed  acoustic  Mach  numbers  less  than 
0.7,  the  overall  sound  pressure  level  radiated  by  an  axisymmetric  jet  is  shown  to  increase  with  the  addition 
of  heat,  and  decrease  above  Ma  >0.7  (Tanna,38  Fisher  et  al.39).  This  is  rather  counter  intuitive  as  the 
dominant  source  term  in  the  cold  jet  pU{Uj,  (Reynolds  stress  term)  would  no  doubt  decrease  in  magnitude 
with  decrease  in  density  p,  resulting  from  the  addition  of  heat.  Ribner,22  after  observing  limited  noise 
reduction  in  jets  of  varying  density  due  to  heating,  suggested  that  the  source  density  played  little  role  in 
affecting  the  overall  noise  of  hot  jets,  and  that  the  entropy  fluctuations  must  dominate  in  these  flows. 

The  addition  of  heat  affects  jets  moving  at  lower  velocity,  through  a  significant  contribution  of  entropy 
fluctuations,  manifest  as  an  increase  in  the  radiated  sound.  However,  at  higher  velocity,  this  effect  is 
diminished  as  the  resulting  overall  sound  pressure  level  is  decreased.  The  Reynolds  stress  term  seemingly 
now  the  driving  mechanism,  as  decreases  in  the  magnitude  of  the  source  density  (in  the  absence  of  additional 
sources),  would  result  in  a  reduction  of  the  radiated  noise.  A  second  interesting  point  of  contention  arises  as 
the  frequency  content  of  the  far-field  spectrum  of  low  Mach  number,  heated  jets  is  examined.  Conversely  to 
that  of  the  cold  jet,  a  larger  angular  distribution  over  which  radiation  occurs  is  seen,  owing  to  an  increased 
magnitude  in  the  refractive  ’cone  of  silence’  region.  The  increase  in  the  overall  sound  pressure  level  is  also 
shown  to  coincide  with  increases  in  the  lower  frequencies  and  a  decrease  in  higher  frequencies.  This  coincides 
with  the  perception  that  the  acoustic  sources  consist  of  quadrupoles  and  dipoles,  with  the  contributions  of 
quadrupolcs  dominant  in  the  cold  jet. 

Fisher  et  al.39  suggested  that  with  the  addition  of  heat,  the  additional  entropy  fluctuation  source  exhib¬ 
ited  a  dominant  contribution  of  dipole-like  sources.  Based  on  these  observations,  Morfey  et  al.40  modeled 
the  the  noise  radiated  from  hot  jets  using  quadrupolc  sources,  with  temperature  fluctuation  terms  of  dipole 
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order,  representing  the  effect  of  heating.  The  resulting  spectrum  showed  differences  in  shape  at  an  angle  of 
90°  with  the  jet  axis,  with  peak  frequency  of  the  dipole  terms  lower  than  that  of  the  quadrupolc.  This  again 
is  consistent  with  the  shift  to  lower  frequency,  previously  exhibited  in  hot  jets,  leading  to  the  fluctuating 
entropy  source  term  (p  -  a^p)  being  dubbed  as  the  ’2nd  dipole  source  term’.  In  Goldstein’s18  formula¬ 
tion  of  the  acoustic  analogy,  it  is  also  noted  that  the  ’enthalpy  flux  term’  corresponds  to  a  dipole  source. 
Viswanathan,41  an  opponent  to  this  view,  studied  extensively  the  affect  of  jet  temperature  on  radiated  noise 
(including  analysis  of  previously  reported  data)  and  concluded  that  the  observed  changes  in  spectral  shape 
were  in  fact  an  effect  due  to  Reynolds  number  and  possible  rig-noise  contamination.  This  contention  among 
researchers  leaves  no  universally  accepted  theory  in  regards  to  the  mechanisms  at  work  in  hot  jets. 


III.  Objective 

The  goal  of  this  effort  has  been  to  identify  the  dominant  source  of  the  radiated  aero-acoustic  noise 
produced  by  high-speed,  heated  jets  and  develop  control  strategies  to  reduce  it.  To  that  end,  two  acoustically 
matched  Mach  0.6  jets,  at  temperature  ratio  Tr=0.93  (cold)  and  Tr— 1.7  (hot)  arc  examined.  The  use  of 
non-intrusivc  Particle  Image  Vclocimctry  (PIV)  to  sample  the  flow,  allowed  a  true  measure  of  the  velocity 
field  to  be  realized  without  fear  of  corrupting  the  radiated  noise  field  intrinsic  to  each  jet.  The  low-order 
modal  dominance  of  each  was  determined  using  Proper  Orthogonal  Decomposition  (POD),  highlighting  the 
low-dimensional  nature  of  this  highly  turbulent  flow  field. 

The  near-field  pressure  signature  surrounding  the  jet  is  sampled  by  an  azimuthal  array  of  fifteen  pressure 
transducers,  at  a  downstream  location  of  2  diameters.  The  predominantly  hydrodynamic  signature  is  shown 
to  correlate  well  with  the  far-ficld  acoustic  response.  The  contribution  to  the  far-ficld  acoustic  spectrum, 
was  determined  by  sampling  the  acoustic  signature  of  each  jet,  using  an  arc  array  of  six  microphones.  These 
quantities  arc  simultaneously  sampled  in  order  to  examine  the  process  by  which  sound  generated  by  the 
axisymmctric  jet  radiates  to  the  far-field.  The  low-dimensional  nature  of  the  source  field  (column  mode  in 
the  hot  jet,  and  helical  in  the  cold),  was  resonated  by  the  near-field  pressure,  and  also  sensed  by  the  far-field 
microphones.  This  consistency  between  all  three  measurements,  source  velocity  field,  near-field  pressure, 
and  far-field  acoustics,  highlights  the  strong  relationship  between  source  generation  mechanism  and  the 
subsequent  propagation  of  fluctuating  pressure,  both  as  hydrodynamic  disturbances  and  far-ficld  acoustic 
disturbances. 

The  Reynolds  stresses  are  weighed  to  glean  some  insight  into  the  entropy  source  term’s  contribution  to 
the  overall  radiated  noise  field.  Where  noise  production  is  concerned,  the  Reynolds  stresses  can  be  viewed  as 
a  simplified  representation  of  the  Lighthill’s  stress  tensor,  Txj  =  puxUj.  Although  not  a  valid  approximation 
to  the  complete  source  term,  with  the  addition  of  heat,  the  weighting  of  these  terms  can  be  telling.  The 
increased  magnitudes  consistently  seen  in  the  hot  jet  would  seem  to  indicate  a  clear  relationship  between 
these  stresses  and  trends  observed  in  overall  sound  pressure  level.  However,  as  this  term  scales  with  density, 
the  vast  difference  in  mean  density  in  each  jet  (pCoid= 1  428  and  phot  =0.759)  demonstrates  that  further 
consideration  of  theses  terms  is  needed.  As  such,  the  Reynolds  stress  terms  contribution  to  noise  generation 
in  the  hot  jet  would  be  diminished  as  compared  to  the  cold  jet,  again  suggesting  the  presence  of  an  additional 
acoustic  source  with  the  addition  of  heat. 

The  acoustic  source  strength  is  sought  in  an  investigation  of  the  quadrupolc  correlations,  following  the 
formulation  of  Proudman68  and  Ribner.69  The  acoustic  power  output,  or  intensity,  in  the  far-ficld  can  be 
expressed  as  fourth  order  space-time  statistics,  calculated  directly  from  the  fluctuating  source  quantities. 
Terms  involving  the  mean  velocity  component  of  the  dominant  flow  direction  in  the  jet  arc  deemed  shear 
noise,  and  those  only  of  fluctuating  quantities,  self  noise.  This  is  more  clearly  articulated  as  noise  contribution 
from  turbulence-mean  flow  interaction  and  turbulence-turbulence  interaction,  respectively.  There  individual 
contribution  is  examined  to  determine  the  dominant  flow  interaction  mechanism  in  each  jet.  As  the  formu¬ 
lation  assumes  isotropic  turbulence  and  neglects  the  contribution  of  the  second  source  term,  (p-a^p),  which 
has  been  shown  to  scale  with  increases  in  temperature  fluctuation,  from  these  results  the  contribution  of  the 
Reynolds  stress  tcrm(s)  is  qualitatively  compared  in  both  jets  to  infer  as  to  the  contribution  of  the  second 
source  term. 
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The  task  of  controlling  the  dominant  sources  of  radiated  noise  is  undertaken,  firstly,  by  modeling  the  flow 
field  dynamics  using  the  full,  unsteady  equations  of  motion  as  a  basis.  PIV  measurements  in  a  longitudinal 
plane  (x,y)  can  be  performed  to  demonstrate  the  evolution  of  flow  structures  as  they  convcct  within  the 
shear  layer.  However,  streamwise  vorticity  is  an  important  quantity  in  the  jet  flow,  being  responsible  for 
much  of  the  dynamics,  and  found  to  be  linked  to  noise  production.17  Therefore,  in  the  present  study 
measurements  arc  performed  in  the  cross-flow  plane  of  the  jet,  i.e.  the  plane  perpendicular  to  the  jet 
centerline.  The  measurement  plane  is  traversed  downstream  by  discrete  increments  to  cover  a  measurement 
volume  extending  from  the  near  jet  exit  (x/D= 3)  to  beyond  the  region  of  collapse  of  the  potential  core 
which  occurs  between  5  and  6  jet  diameters  downstream.  At  this  location,  turbulence  intensity  reaches  a 
maximum  and  the  near-field  pressure  dynamics  seem  to  relate  best  with  the  far-ficld  sound,  as  shown  by 
Hall  et  a/.52  .  To  capture  and  rebuild  the  dynamics  of  the  flow,  a  low-dimensional  experimental  approach  is 
taken  here,  in  which  a  large  database  of  the  acceleration  field  (and  thereby  the  velocity  field)  is  needed.  Since 
the  velocity  data  contains  the  information  described  by  the  Navicr-Stokcs  equations,  a  system  of  ordinary 
differential  equations  (ODE’s)  is  fit  to  the  data  to  provide  a  dynamical  system  capable  of  predicting  the 
low  order  dynamics  of  the  jet  flow  field,  given  an  initial  condition.  This  capability  carries  interest  for  flow 
control  applications  where  a  dynamical  estimate  of  the  state  of  the  flow  can  be  of  great  use  in  conjunction 
with  direct  measurement. 


IV.  Low-dimensional  Jet  Analysis 

In  many  aspects,  research  involving  the  complex  flows  associated  with  turbulence  has  been  limited  by 
the  development  of  measurcmcnt/computational  tools.  The  complexity  of  the  full  non-linear  equations  of 
motion  makes  numerical  results  such  as  direct  numerical  simulations  (DNS)  limited,  even  with  improved 
processing  power  and  speed  of  computers.  Experiments,  relying  on  measurement  prove  futile  in  many  cases, 
as  necessary  flow  properties  arc  at  times  impossible  to  acquire  without  corruption.  Leaving  techniques  such 
as  those  incorporating  modeling  and  statistics  at  the  forefront  of  the  development  of  new  theories.  These  set 
of  experiments  couples  measurements  and  statistics  to  glean  insight  into  the  source  velocity  characteristics, 
near-field  pressure  signature,  and  acoustic  features,  distinct  to  both  the  cold  and  hot  jet. 


A.  Experiment 

Two  separate  sets  of  experiments  were  carried  out  at  Mach  0.6  exit  velocity,  with  varying  core  flow  tempera¬ 
tures,  to  quantify  the  effect  of  elevated  temperature  on  the  acroacoustic  sources.  The  first  at  a  temperature 
of  20 °C  (  293 °K)  which  will  be  referred  to  as  the  cold  jet  and  a  second  at  243°C  (  516 °K)  referenced  as 
the  hot  jet.  In  both  studies,  the  ambient  conditions  were  held  constant  at  T=20°C  and  p=16ar. 


1 .  Facility 

All  experiments  are  conducted  in  the  Syracuse  University  Skytop  Ancchoic  Chamber,  located  at  Syracuse 
Univcrsitys  Skytop  campus.  The  facility,  is  approximately  30  years  old,  having  been  constructed  in  the  early 
seventies  under  the  leadership  of  Dosanjh  et  al47  and  designed  specifically  with  the  purpose  of  conducting 
acroacoustic  studies.  All  exterior  walls,  flooring  and  ceiling  arc  constructed  of  reinforced  30.5cm  thick 
concrete.  The  chamber,  encompassing  206m3,  is  acoustically  treated  using  high-density  fibreglass  wedges  on 
all  interior  walls,  with  a  cut-off  frequency  of  150 Hz.  The  interior  of  the  chamber  can  be  seen  in  Figure  1.  A 
detailed  description  of  the  facility  and  it’s  outfitted  equipment  is  presented  by  Tinncy  et  al 48 

The  experimental  jet  rig  in  Figure  2  is  comprised  of  three  main  parts,  two  stainless  steel  sections,  joined 
at  a  flange  connection,  and  a  machined  nozzle.  The  nozzle  interior  utilizes  a  matched  5th  order  polynomial 
profile  for  the  contraction  from  15.2cm  pipe  to  a  5.08cm  exit  diameter.  The  mounting  stand,  fashioned  from 
carbon  steel  I-beams  and  designed  to  allow  for  expansion  of  piping  at  elevated  temperatures,  completes  the 
experimental  jet  rig.  The  stand  compensates  for  this  growth  using  several  high  strength  linear  rails  providing 
a  rigid  single  degree  of  freedom  system,  in  the  direction  of  flow. 
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Figure  1.  Interior  of  Syracuse  University  Skytop  Anechoic  Chamber. 


The  power  plant  of  the  facility  is  comprised  of  a  100/ip,  two-stage  (2  piston)  Joy  compressor,  housed  in 
a  separate  building.  Operating  as  a  blow-down  facility,  the  compressed  air  is  stored  in  five  tanks  with  a 
volumetric  capacity  of  6.2 m3  each,  to  a  maximum  pressure  of  approximately  3400 kPa.  Once  at  pressure, 
desired  exit  flow  rates  of  the  jet  rig  can  be  achieved.  This  is  done  through  the  use  of  a  robust,  industrial 
grade  Allen  Bradley ,  Pro-Logix  Control  (PLC)  system,  which  drives  the  flow  process  and  monitors/controls 
much  of  the  equipment  ( i.e .  industrial  make-up  air  unit,  electric  heater,  flow-sense  switch,  numerous  pressure 
transducers,  compressed  air  tank  pressure  and  throttle,  and  exhaust  fan).  As  input  sources  to  the  jet  rig 
operation  arc  the  storage  tanks’  pressure,  static  pressure  in  the  nozzle,  and  both  plenum  and  chamber, 
pressure  and  temperature.  From  a  dedicated  PC,  the  supply  tanks  pressure  can  be  monitored  until  sufficient 
pressure  is  reached.  Once  at  pressure,  a  user  defined  exit  velocity,  in  terms  of  Mach  number,  is  achieved 
based  on  the  ratio  of  static*  pressure  in  the  nozzle  to  chamber  pressure.  The  pneumatically  actuated  valve 
controlling  the  air  supply  is  then  throttled  to  a  corresponding  percentage  opening,  giving  the  prescribed  exit 
velocity.  The  operational  time  of  the  jet  facility,  based  on  maximum  pressure,  is  just  below  twenty  minutes 
when  discharging  at  Mach  0.85,  for  the  50.8mm  nozzle. 


Governing  the  flow  process,  and  maintaining  conditions  during  experiments,  arc  a  scries  of  inputs  from 
pressure  transducers  and  thermocouples,  used  to  monitor  the  flow’s  properties.  The  pressure  ratio  used 
to  control  flow  rate  is  given  by  a  set  of  in-pipe  and  chamber  transducers.  This  ratio  drives  a  diaphragm- 
actuated,  sliding-stem  control  valve  to  maintain  a  desired  exit  velocity.  The  static  pressure*  used  in  the  ratio 
to  determine  the  exit  Mach  number,  has  been  shown  by  Tinney  et  a/.,48  to  be  a  valid  approximation  to  the 
total  pressure  used  in  the  iscntropic  relationship, 


Po 

P 


7-1 

2 


7/7-1 


(10) 


where  p0  and  p  are  total  and  ambient  pressure,  respectively.  The  calibration  between  the  nozzle  static 
pressure  (transducer)  and  total  pressure  (pitot  probe)  reveals  a  nearly  linear  relationship  between  the  two, 
exhibiting  minimal  differences.  The  ratio  of  pipe  static  pressure  to  ambient  can  therefore  be  used  in  place 
of  total  to  ambient  as  specified  by  Eq.  10  above. 


2.  Instrumentation 

The  near-field  pressure  is  acquired  using  an  azimuthal  array  of  fifteen  Kulite  model  XCE-093-5G  transducers. 
These  arc  positioned  in  a  array  at  two  diameters  downstream  from  the  jet  exit,  equally  spaced  at  24°  apart, 
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Figure  2.  Experimental  jet  rig  (left)  left-side  view  (right)  right-side  schematic  drawing. 


and  approximately  1cm  outside  the  spreading  shear  layer.  The  array  is  mounted  to  the  nozzle  using  a  cap 
fitted  with  a  scries  of  holders.  Their  orientation  is  pictured  in  Figure  3.  The  frequency  response  of  these 
transducers  is  150 kHz,  with  dynamic  range  DC  to  50kHz,  and  a  20mV /6.8kPa  sensitivity.  Their  small 
size,  (shaft  diameter  2.36xl0_2cm)  makes  them  an  ideal  measurement  tool,  evoking  minimal  obstruction  of 
flow  characteristics.  These  are  continuously  sampled  at  40,960 Hz  and  sub-sampled  during  analysis  in  blocks 
of  8192  samples,  giving  a  frequency  resolution  of  df=bHz. 


Figure  3.  Azimuthal  array  of  Kulite  transducers. 


To  sample  the  far-ficld  acoustic  characteristics,  six  G.R.A.S.  Type  40 BE  prepolarized,  frcc-field  micro¬ 
phones  are  used.  The  frequency  response  is  4 Hz  to  100 kHz(-\-/  -  1  dB)  or  10 Hz  to  40kHz(+/  -  2dB),  with 
dynamic  range  \66dB(ref20pPa)  and  a  sensitivity  of  4 mV/ Pa.  These  front-vented  microphones  are  capa¬ 
bly  of  operating  at  temperatures  ranging  from  -5 °C  to  120°C.  Excitation  is  provided  by  a  G.R.A.S.  Type 
26CB  preamplifier  with  a  frequency  response  range  between  4 Hz  to  100kHz  with  an  uncertainty  of  ±0.2 dB. 
These  arc  positioned  on  a  boom  marking  an  arc  fixed  at  85  jet  diameters  from  jet  centerline.  The  array  of 
microphones  are  then  mounted  on  threaded  rods  extending  from  the  boom,  placing  there  diaphragms  on  an 
arc  of  75  diameters.  The  array  is  pictured  in  Figure  4,  with  its  orientation  also  diagramed.  The  downstream 
microphone  (mic.  6)  is  positioned  at  an  angle  of  0=15°  from  the  jet  centerline  and  arrayed  at  15°  thereafter, 
with  the  sixth  (mic.  1)  at  0=90°  from  jet  exit. 

Heating  the  jet  core  flow  is  achieved  using  a  470kW  Chromalox  electric  circulation  heater.  The  heater 
utilizes  of  scries  of  sheathed  rods  that  power  on  and  off,  as  needed,  to  achieve  the  desired  user  input 
temperature.  It’s  operation  is  controlled  by  the  PLC  system,  with  a  necessary  requirement  being  the  output 
of  a  flow-sense  switch,  to  prevent  damage  to  the  rods.  The  output  temperature,  measured  by  a  thermocouple 
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Figure  4.  Boom  array  of  far-field  microphones. 


upstream  of  the  jet  rig,  is  a  function  of  a  number  of  components.  The  outside  air  temperature  that  supplies 
the  compressor,  flow’s  properties  as  it  exits  the  storage  tanks,  and  the  initial  temperature  of  the  heater  all 
play  a  role.  As  the  flow  is  heated  through  convection,  this  delicate  balance  inhibits  a  user  prescribed  output 
temperature,  but  rather  only  allows  a  heater  coil  temperature  to  be  defined.  A  nozzle  exit  flow  temperature 
of  53 8°C  (1000°F)  can  be  achieved,  at  a  flow  rate  of  41.3m3/mm. 


Figure  5.  FIV  setup,  showing  camera  orientation  and  laser  positioning. 


A  Dantec  Dynamics  3-componcnt,  stereoscopic  PIV  system  is  used  to  acquire  the  jet’s  velocity  field.  The 
system  incorporates  two  HiSense  8-bit  resolution  CCD  cameras  (1280x1024  pixels),  a  pair  of  pulsed  New 
Wave  Research  200 mJ  Nd:YAG  lasers,  and  an  optical  laser  sheet  generator.  The  cameras  arc  positioned  at 
an  angle  of  45°  with  jet  centerline,  mounted  just  upstream  of  the  jet  rig  (at  a  height  inline  with  the  piping), 
on  cither  side.  The  laser  sheet  illuminates  the  field  from  beneath  in  the  radial-azimuthal  (r  -  6)  plane, 
perpendicular  to  the  jet  exit.  This  orientation,  shown  in  Figure  5,  establishes  the  streamwise  direction  (x), 
as  the  out-of-planc  component  of  the  velocity  field.  The  iterative  correlation  scheme  used  to  resolve  the 
image  pairs  captured  by  each  camera,  gives  in  plane  particle  displacement,  and  is  thus  a  two  component 
representation  of  the  spreading  shear  layer.  The  single  resultant  velocity  maps  (one  from  each  camera)  arc 
transformed  to  realize  the  third  component,  based  on  out-of-planc  displacement.  This  orientation  allows  a 
greater  limit  of  measurable  velocity  range,  as  the  time  step  between  laser  pulses,  is  now  not  a  direct  function 
of  jet  streamwise  velocity.  Acquiring  the  strong  mean  flow  in  this  direction  would  directly  limit  the  system, 
as  particle  displacement  within  a  given  interrogation  area  would  be  large;  requiring  very  short  duration 
between  captured  images.  In  this  high  speed,  Mach  0.6  flow,  a  duration  of  5//s  between  pulses  was  used 
resulting  in  particle  displacement  on  the  order  of  a  millimeter  (laser  sheet  thickness  of  5mm).  The  system’s 
positioning  is  fixed,  including  both  cameras  and  laser  sheet  generator,  and  mounted  on  a  single  degree  of 
freedom  traverse.  As  an  initial  calibration  locates  camera  positioning  with  reference  to  laser  sheet,  this  limits 
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the  need  to  recalibrate  the  system  at  different  downstream  sampling  locations. 

A  key  component  in  studies  involving  the  use  of  optics  based  measurement  tools,  is  seeding.  To  seed 
the  core  flow,  a  PIVTEC,  twelve  Laskin  nozzle  seeder,  with  olive  oil  as  its  fuel,  is  used  in  the  cold  study. 
Its  design  allows  pressures  upwards  of  lObar,  well  above  the  back  pressure  experienced  from  injecting  into 
the  jet  piping.  Median  particle  sizing  is  on  the  order  of  1/im.  For  the  hot  study,  a  seeder  was  designed 
in-house,  using  a  3M  ceramic  microsphere  particulate  seeding  medium.  The  pressure  vessel  included  two 
angled  injection  ports  to  promote  mixing  and  a  fine  10 fim  mesh  before  the  exit  port  to  prevent  particles 
from  clumping.  An  ad  hoc  investigation  of  the  mean  particle  sizing  revealed  a  value  on  the  order  of  1-10 fim. 
The  coflow  is  also  seeded  in  both  studies,  providing  a  uniform  distribution  of  particles  in  the  entrainment 
region  as  well.  A  TSI  model-9307  oil  droplet  generator  is  used,  first  filling  the  plenum,  then  pushed  by  the 
make-up  air  unit  into  the  circular  co-flow  duct  surrounding  the  jet  piping.  This  seeder  also  uses  olive  oil 
as  its  seeding  medium.  Based  on  manufacture  specifications,  the  median  particle  size  on  the  order  of  1  fim. 
These  sizes  are  optimal  based  on  discussions  of  Mclling49  on  the  effect  of  tracer  particle  sizing  and  other 
aspects  relevant  to  optical  measurements  tools. 

The  sampling  of  the  instrumentation  mentioned  above  is  accomplished  through  the  use  of  a  National 
Instruments  PXI  system,  with  a  866MHz  NI-8175  embedded  controller.  The  system  houses  three  NI-4472 
boards,  each  with  eight  available  24-bit  resolution  channels,  and  all  with  independent  analog  to  digital 
converter.  Of  the  twenty- four  available,  fifteen  channels  are  dedicated  to  the  Kulite  transducers,  six  to  the 
G.R.A.S ,  and  one  to  the  timing  of  the  Dantec  Dynamics  PIV  system.  The  system  is  capable  of  continuous 
sampling,  streaming  all  channels  simultaneously  at  rates  above  50 kHz.  The  sampling  rate  of  40,960 Hz  used 
falls  well  within  these  limits.  A  Krohn-Hite  model  filter  is  used  as  a  signal  conditioner  before  acquisition 
of  the  far- field  microphones.  The  filters  attenuate  at  a  rate  of  24 dB/ octave  and  have  a  maximally  flat 
(Butterworth)  amplitude  response.  The  high  pass  setting  of  100 Hz,  is  intended  to  eliminate  possible  low- 
frequency  vibrations  from  the  boom  array.  A  20 dB  gain  is  also  applied. 


B.  Cold  study 


In  an  attempt  to  examine  the  effect  of  temperature  on  the  acroacoustic  sources  of  the  jet,  a  cold  study  was 
first  performed  to  serve  as  a  ’control’  case.  At  an  exit  Mach  number  of  0.6,  the  core  flow  was  held  constant 
at  approximately  20°C.  Based  on  isentropic  flow  relationships,50 
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the  nozzle  total  to  static  pressure  ratio  (NPR)  is  1.2755,  with  temperature  ratio  of  1.072.  However,  as  is  the 
ease  in  jet  acoustic  studies,  a  new  temperature  ratio  and  Mach  number  arc  defined  and  will  be  referred  to  as 
the  acoustic  temperature  ratio  (Tr)  and  Mach  number  (Ma).  Where  Tr  is  defined  by  the  ratio  of  jet  total 
temperature  at  exit  over  the  ambient  temperature  and  Ma  as  the  ratio  of  jet  exit  velocity  to  ambient  sound 
speed  (cioo).  For  completeness,  the  jet  Mach  number  is  expressed  as  My 
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indices  referenced  as  O-(total),  j-(jet),  oo-(ambient) 

The  relevant  properties  of  the  cold  study  arc  itemized  in  Table  4  below.  Velocity  measurements  arc  acquired 
at  half  diameter  increments  from  4.5-10D,  with  ambient  conditions  held  constant  at  T=20°C  and  p=16ar, 
by  use  of  the  co-flow  system.  Throughout  the  remainder  of  the  document,  this  study  is  referenced  as  the 
cold  jet. 
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Table  1.  Cold  Study 


Too 

Tr 

Ma 

M, 

NPR 

Re 

20 °C  (  293 °K) 

0.93 

0.6 

0.6 

1.2755 

6.9xl05 

C.  Hot  study 

In  the  second  set  of  experiments,  the  jet’s  core  flow  temperature  is  increased  to  243 °C  (  516 °K).  As  the 
acoustic  properties  in  the  jet  are  a  function  of  both  exit  velocity  and  temperature,  the  exit  velocity  (Uj)  must 
be  held  constant  to  isolate  the  effect  of  temperature.  The  properties  of  air  change  at  elevated  temperature, 
so  to  maintain  a  constant  exit  velocity  translates  into  lowering  the  jet  Mach  number;  whereby  the  sound 
speed  in  Eq.15  increases  with  increases  in  temperature,  as  a  result  of  a  decrease  in  fluid  density.  In  this 
instance,  it  is  prudent  to  use  the  acoustic  Mach  number,  coupled  with  a  temperature  ratio  as  a  measure  of 
the  jet’s  exit  conditions.  The  study  is  thus  described  as  an  investigation  of  the  acroacoustic  sources  of  an 
’acoustically  matched’  Mach  0.6  jet,  at  varying  temperature  ratio. 

The  air  exiting  the  compressor,  stored  in  the  large  tanks  powering  the  jet,  must  be  heated  upstream 
of  it’s  exit.  As  it  convects  over  the  sheathed  rod  of  the  electric  heater,  only  a  coil  temperature  can  be 
determined.  Based  on  the  velocity  of  the  air  in  the  upstream  piping,  a  setting  of  approximately  315°C  was 
necessary  to  achieve  an  average  temperature  of  243°C,  over  the  duration  of  a  single  run  of  approximately  16 
minutes.  In-pipe  temperature  is  measured  by  the  thermocouple  positioned  just  upstream  of  a  honeycomb 
flow  straightener,  housed  within  the  connection  of  piping  and  nozzle.  The  relevant  properties  of  the  hot  study 
arc  itemized  in  Table  2  below.  The  jet  Mach  number  now  falls  to  0.44  and  Reynolds  number  almost  half  of 

Table  2.  Hot  Study 


Too 

Tr 

Ma 

Mj 

NPR 

Re 

243°C  (  516° A") 

1.7 

0.6 

0.44 

1.142 

2.8xl05 

the  cold  study,  however  remaining  a  turbulent  flow  field  still  within  the  compressible  flow  regime.  Velocity 
measurements  arc  again  acquired  at  half  diameter  increments  from  4.5-10 D,  with  ambient  conditions  held 
constant  at  T=20°C  and  p=lbar,  by  use  of  the  co-flow  system.  This  study  is  hereafter  referenced  as  the 
hot  jet.  Any  increase  in  chamber  temperature  or  pressure  was  indiscernible  to  measurement  tools,  i.e. 
thermocouple  positioned  at  the  base  of  the  jet  rig  and  transducer  on  chamber  floor.  There  was  also  no 
noticeable  expansion  of  the  jet  piping  under  the  described  thermal  loading. 


D.  Near-field  pressure 

The  averaged  power  spectral  density  of  the  sampled  fluctuating  pressure  surrounding  the  cold  jet  is  now 
examined. 

s^u)  =  f  (pj(/)p;(/))  (i6) 


(*)  denotes  the  complex  conjugate 

Similar  shapes  arc  echoed  by  each  transducer  as  a  transfer  function,  derived  based  on  the  averaged  amplitude 
spectrum  over  all  transducers,  is  applied  prior  to  the  following  analysis.  The  similar  orientation  of  the  trans¬ 
ducers  and  the  mean  azimuthal  invariance  of  the  axisymmetric  jet,  dictates  that  all  averaged  power  spectral 
densities  will  be  identical.  The  application  of  the  transfer  function  has  negligible  impact,  as  the  transformed 
signal  simply  experiences  an  amplitude  modulation,  retaining  all  phasing  and  frequency  characteristics. 

The  data  is  analyzed  using  methodology  presented  in  the  Arndt  et  al.29  investigation  of  the  pressure  field 
surrounding  the  turbulent  jet.  A  simple  model  based  on  a  point-source  solution  to  the  spherical  wave  equation 
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was  utilized  in  order  to  divide  the  pressure  fluctuations  into  two  forms:  propagating  or  acoustic  fluctuations 
(in  phase  with  the  velocity  fluctuations),  and  non-propagating  or  hydrodynamic  fluctuations  (out  of  phase 
by  90°  with  the  velocity).  Using  the  assumption  that  long  wavelength  disturbances  arc  associated  with  large 
sources  and  that  short  wavelength  disturbances  are  associated  with  small  sources,  an  expression  is  derived 
for  the  intensity  of  the  hydrodynamic  field  (for  a  constant  wave  number  u  =  a0k )  which  would  have  a  spatial 
decay  of  /  oc  r-6,  where  r  denotes  a  radial  distance  outward  from  the  center  of  the  source  (mixing  layer 
centerline).  While  at  constant  r,  in  the  inertial  subrange,  the  spectral  variation  would  be  something  of  the 
order  of  /  oc  k~2^3k~6  oc  k~6  67.  More  importantly,  it  is  demonstrated  that  the  transition  between  these 
two  slopes  (the  hydrodynamic  near  field  and  the  acoustic  far  field  regions  of  the  pressure)  was  frequency 
dependent  and  occurred  predominantly  around  kr= 2,  exhibiting  a  /  oc  k~ 2  slope.  This  analysis  has  been 
revisited  more  recently  by  Hall  et  al 52  and  Coiffet  et  a/.51  who  revealed  strong  interaction  mechanisms 
between  the  reactive  and  propagative  components  of  the  pressure  field  just  outside  of  the  rotational  mixing 
layer  region  of  a  Mach  0.85  &  0.30  round  jet,  respectively. 


Figure  6.  Transformed  near-field  pressure  spectra,  x/D=  2,  r/D  ~1  ( all  15  transducers  lay  ontop  each  other). 


The  non-dimensional  spectral  densities,  normalized  by  dynamic  head  (O.bpU2)2,  arc  presented  in  Figure  6. 
These  arc  averaged  over  200  blocks  of  8192  data  points,  equivalent  to  40  seconds  of  sampled  data,  at  an 
acquisition  rate  of  40,960 Hz.  Unless  otherwise  noted,  all  presented  spectra  are  analyzed  using  these  criteria. 
The  slope  exhibits  a  roll-off  consistent  with  that  of  hydrodynamic  spectral  variations,  and  does  manifest  a 
bifurcation  in  concurrence  with  Arndt  et  al.29  however  at  kr  =  3,  not  kr  =  2.  This  has  been  the  matter  of 
some  contention  among  researchers.  In  the  plot  presented,  there  is  also  only  a  small  region,  following  the 
kr  =  3  separation  which  exhibits  a  slope  near  kr~2.  The  positioning  of  the  transducers  ( x/D—2  and  r/D  ~  1, 
just  outside  the  spreading  shear  layer)  may  play  a  role  here,  as  Ko  &:  Davies30  and  Lau  et  al.5  have  reported 
evidence  regarding  the  absence  of  visible  vortices  or  the  radially  inward  entrainment  of  velocity  in  the  near 
lip  region  of  the  flow  field.  Suggesting  that  the  development  of  the  source  field  at  this  location  is  premature. 
In  separate  measurements  acquired  at  a  downstream  location  x/D=S ,  and  several  radial  locations  (r/D=1.9 
to  2.68),  this  dissemblance  between  the  reactive  (hydrodynamic)  and  propagative  (acoustic)  components 
of  the  pressure  field  is  more  clearly  articulated.  This  location  is  known  to  lie  within  the  dominant  sound 
production  region,  exhibiting  a  larger  amount  of  energy  with  k~ 2  roll-off  slope,  indicating  the  presence  of  an 
acoustic  source.  This  analysis  dictates  that  care  must  be  taken,  where  sensing  of  the  deterministic  sources 
of  sound  near  the  jet  exit  is  concerned.  At  the  present  measurement  location  of  x/D= 2,  although  seemingly 
unlikely  to  apprehend  a  realizable  acoustic  signature  until  some  distance  downstream  where  the  propagating 
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acoustic  signature  is  dominant,  x/D=6-10,  it  is  demonstrated  that  much  information  can  be  gleaned  from 
this  near  lip  region.  This  is  examined  in  greater  detail  in  the  coupling  of  the  near  and  far-ficld  pressure. 

Also  clearly  illustrated  in  Figure  6,  is  the  marked  increase  in  low  frequency  energy  seen  in  the  fluctuating 
pressure  spectra  of  the  hot  jet,  accompanied  by  a  slight  decrease  at  high  frequency.  With  the  lowered  Re 
of  the  hot  jet,  2.8xl05  from  6.9xl05,  it  may  be  expected  that  less  fluctuations  would  exist  as  the  jet  moves 
towards  more  laminar  behavior.  However,  an  increase  in  the  pressure  fluctuations  is  observed,  prms  =0.3163 
in  the  cold  jet  and  prms=0.4179  in  the  hot.  It  is  well  known  that  due  to  increased  spatial  growth  rate  of 
instability  waves  in  hot  jets,  increased  mixing  is  manifest.  This  feature  is  reflected  here  in  the  spectra,  by  an 
increase  in  energy.  The  phenomenon  observed  suggests  that  there  is  an  increased  low  frequency  contribution 
of  the  dominant  structures  near  the  jet  exit  in  the  hot  jet,  which  is  sensed  by  the  surrounding  pressure 
near-field. 

Additional  features  of  the  pressure  field  are  shown  in  Figure  7,  comprising  the  coherency  spectra  between 
the  azimuthal  array  of  probes  of  both  the  cold  and  hot  jets.  The  spectra  is  averaged  over  all  azimuthal 
separation,  as  correlations  arc  evaluated  with  each  probe  as  the  0  =  0  coordinate  and  then  averaged  over 
similar  AO  separation.  The  plots  only  show  the  first  180  degrees  of  separation,  as  the  coherency  at  greater 
AO  separation,  lie  on  similar  curves.  The  curves  peak  at  small  azimuthal  separation  (AO  =  24),  then  decrease 
with  increasing  separation,  suggesting  the  presence  of  a  dominant  Fourier-azimuthal  helical  mode  structure. 
Similar  coherency  is  exhibited  at  nearly  78%  in  the  cold  jet,  and  just  over  72%  in  the  hot  jet.  However,  at 
larger  AO  separation,  the  hot  jet  exhibits  a  stronger  correlation,  more  indicative  of  an  axisymmctric  column 
mode  structure.  This  trend  suggests  the  presence  of  differing  dominant  modal  contribution  in  each  jet. 
There  is  also  a  marked  difference  in  the  dominant  frequency  at  which  these  peaks  occurs.  In  the  cold  jet 
(7a),  a  more  broadband  peak  persists,  centered  at  approximately  Stp= 0.5.  The  hot  jet  (7b),  exhibits  a  more 
distinct  peak,  over  all  separation,  at  a  lower  Stp  near  0.3,  suggesting  that  a  greater  distribution  of  scales  is 
present  in  the  cold  jet. 

To  further  investigate  the  observed  trends,  the  distribution  of  the  modal  characteristics  of  the  near-field 
pressure  correlation  is  presented.  The  cosine  transform  is  not  utilized,  for  purposes  of  computational  ease 
and  time.  Instead,  the  real-part  of  the  cross-spectra,  Fourier  transformed  in  azimuth  (AO  — ►  m),  is  presented 
( only  positive  modes,  +m).  Striking  in  Figure  8,  is  the  scarcity  of  higher  azimuthal  modal  energy  that  has 
been  shown  to  govern  the  primary  features  of  the  turbulent  source  field,  as  only  the  first  two  modes  (0&1) 
exhibit  a  significant  amount  of  energy,  in  both  the  cold  and  hot  jet.  This  modal  behavior  has  also  been 
shown  by  Petersen,9  Arndt,29  and  Hall  et  a/.54 

The  azimuthal  mode  spectra  illustrates  a  similar  picture  as  did  the  coherency,  where  the  dominant 
peak  frequencies  in  the  cold  jet  span  a  greater  range  than  seen  in  the  hot  jet.  Mode  0,  appears  to  peak 
at  approximately  St £>=0.4,  then  just  below  that  value  in  the  hot  jet.  This  trend  is  echoed  at  higher 
mode  number.  In  the  modal  distribution  of  both  jets,  the  column  mode  peaks  at  higher  frequency,  then 
a  progressively  lower  frequency  is  exhibited  at  each  consecutive  mode  number.  Although  more  broadband 
than  its  hot  counterpart,  the  peak  observed  in  the  mode  0  spectra  exhibits  a  more  distinct  shape,  when 
compared  to  the  broad  peaks  exhibited  at  higher  mode  number,  in  both  jets.  It  is  of  note  to  highlight  that 
the  energy  content  of  each  docs  however  differ  greatly,  with  mode  0  of  the  hot  jet  more  than  doubling  that 
of  the  cold. 

The  total  azimuthal  modal  energy  content  is  presented  in  Figure  9,  where 

Fpp  =  ^  Spp(m,  f)df  (17) 

Conclusions  can  now  be  drawn  as  to  the  distinct  dominant  structures  present  in  each  jet.  The  energy 
distribution  is  summed  over  positive  and  negative  modes  (±m)  here,  as  to  only  view  the  positive  contribution 
of  the  modes  would  give  misleading  information.  The  axisymmctric  column  mode  is  shown  to  dominate  the 
flow  field  in  Figure  8,  however  it  is  demonstrated  here,  as  suggested  by  the  coherency  plots,  that  the  second 
mode  (helical  mode)  has  a  significant  contribution  as  well.  This  is  especially  noticeable  in  the  case  of  the 
cold  jet  Figure  9a,  where  the  mode  0  &  1  contributions  arc  nearly  identical.  In  the  hot  jet  there  is  a  marked 
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distinction  between  the  contribution  of  these  two  modes.  Modes  0  &:  1,  account  for  nearly  75%  of  the  total 
energy  in  both  the  cold  and  hot  jet,  with  the  hot  jet  exhibiting  40%  more  overall  energy.  The  distribution 
in  each  case  is  highlighted  in  Table  3  below,  suggesting  that  the  presence  of  the  stronger  helical  mode  in 
the  cold  study,  may  gives  rise  to  the  broadening  of  peaks  exhibited  in  the  coherency  spectra  (Figure  7a). 
The  more  even  distribution  over  the  modes  suggests  a  greater  distribution  in  the  dominant  scales.  As  such, 
the  narrow-band  peak  present  in  the  coherency  curves  of  the  hot  jet,  coinciding  with  a  dominant  mode  0 
contribution,  dictates  the  presence  of  a  consistently  larger  sealed  structure. 


Table  3.  Modal  energy  distribution  ( percentage  of  total). 


Cold 

Hot 

mode  0:  38% 

mode  0:  45% 

mode  1:  37% 

mode  1:  31% 

From  this,  it  can  also  be  inferred  that,  it  is  the  presence  of  the  helical  mode  which  significantly  diminishes 
the  strength  of  the  coherency  at  larger  separation  (A0=48-168).  The  hot  jet  exhibits  greater  coherency  at 
these  separation,  as  the  column  mode’s  contribution  is  more  pronounced  here  (an  entirely  helical  structure 
would  be  largely  uncorrelated  at  180  degrees  of  separation,  where  as  a  column  structure  would  be  well  correlated 
over  all  separation).  The  low  frequency  energy  increase  seen  in  the  fluctuating  pressure  power  spectral 
density  can  therefore  be  linked  to  the  presence  of  a  dominant  axisymmctric  column  mode  0.  The  peak  value 
of  Sto=0A  seen  in  the  mode  spectra  is  consistent  with  the  low  frequency  peak  observed  in  both  the  power 
and  coherency  spectra  presented. 


E.  Far-field  acoustic  signature 

Similarly  to  the  investigation  of  the  near-field  pressure,  the  spectral  densities  of  the  far-field,  or  acoustic, 
pressure  signature  arc  presented  in  Figure  10,  though  as  a  function  of  wavenumber,  not  kr.  (//mm(/j  is 
calculated  in  similar  fashion  as  Spp(f)  using  Eq.  16).  The  radial  location  is  neglected  as  only  the  character¬ 
istics  of  the  spectral  variations  in  the  signal  arc  of  concern.  As  can  be  expected,  the  curves  depict  a  slope 
primarily  scaling  with  a  slower  decay  of  /  a  k~2 ,  consistent  with  that  of  a  dominant  acoustic  field. 

In  both  the  cold  and  hot  jet,  at  shallow  aft  angles  (15-30  degrees  from  jet  centerline,  corresponding  to 
microphones  5  &:  6)  the  signal  exhibits  larger  magnitude,  than  at  large  aft  angles  (75-90  degrees  corresponding 
to  microphones  1  &:  2).  This  is  consistent  with  the  directional  nature  of  the  jet’s  far-field  radiated  spectrum. 
As  the  noise  source  is  deemed  as  the  interaction  of  turbulent  fluid  convecting  with  the  bulk  mean  flow,  the 
distinct  far-field  acoustic  spectrum  of  the  jet  is  a  result  of  a  distorted  radiated  field  from  a  point  source.  The 
magnitude  of  the  resulting  field  is  enhanced  in  the  direction  of  convection.  Where  the  hot  jet  is  concerned, 
there  is  a  noticeably  larger  increase  exhibited  over  the  range  of  microphones,  seen  predominantly  at  shallow 
aft  angles,  with  only  the  downstream  microphones  positioned  from  0=15-45°,  exhibiting  significant  change. 
As  jet  noise  is  known  to  scale  with  velocity,  at  a  fixed  exit  velocity  ( MaCold ^  =  A faHot^)  between  the  two 
experiments,  this  increase  can  be  attributed  to  a  marked  contribution  of  acoustic  energy  from  an  additional 
( entropy )  source  due  to  heating,  in-line  with  the  theory  of  hot  jets.  This  is  in  direct  relation  to  the  diminishing 
of  the  Reynolds  stress  source  term  (pUiUj)  with  decreases  in  fluid  density.  The  normalization  by  (0 .bpU2)2 
carried  out  here  is  similar  to  that  of  Mollo-Christcnsen  &  Narasimha,55  who  used  the  jet’s  dynamic  head  to 
establish  a  non-dimensional  similarity  relationship  for  the  spectra  of  emitted  sound.  As  the  Reynolds  stress 
source  term  can  be  thought  to  scale  with  (pU2),  the  normalized  spectra  presented  can  ultimately  be  thought 
of  as  depicting  only  the  contribution  of  the  second  entropy  source  term  (p  -  a%p). 

In  this  regard,  the  spectra  presented  in  Figure  11a  highlights  this  contribution  in  the  cold  and  hot  jet  at 
microphone  locations  0=15°  &:  90°.  There  is  a  clear  increased  contribution  in  the  hot  jet  at  low  frequency, 
and  a  decrease  at  high  frequency,  with  a  maximum  difference  seen  at  Stp= 0.2,  at  the  0=15°  microphone 
location.  This  is  consistent  with  the  low  frequency  energy  of  the  dipole  like  distribution  of  the  second 
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(b) 

Figure  7.  Coherency  spectra  over  first  180  degrees  of  separation.  (a)Cold  jet.  (b)Hot  jet. 
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(a) 


(b) 

Figure  8.  Azimuthal  mode  contribution  to  cross-spectra  ( modes  0-3).  (a)Cold  jet.  (b)Hot  jet. 
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(a) 


(b) 


Figure  9.  Modal  energy  content  ( summed  over  ±modes).  (a)Cold  jet.  (b)Hot  jet. 
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Figure  10.  Normalized  far-field  spectra.  (a)Cold  jet.  (b)Hot  jet. 
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source  term  as  modeled  by  Morfey  et  al40  The  increase  as  a  measure  of  sound  pressure  level  is  depicted  in 
Figure  lib.  There  is  a  greater  disparity  in  noise  level  seen  between  the  two  microphones  in  the  hot  jet,  as 
the  cold  jet  exhibits  only  a  change  of  ~10 dB  and  hot  jet  ~14 dB,  over  the  range  of  aft  angles  from  0=15° 
through  0=90°.  Again,  the  maximum  change  occurs  near  the  Sto= 0.2  value. 

The  decrease  in  high  frequency  contribution  seen  in  the  comparison  of  the  two  microphone  locations 
(0=15°  k  0=90°),  exhibiting  the  maximum  and  minimum  change,  respectively,  suggests  that  the  apparent 
increase  in  acoustic  energy  may  more  correctly  be  described  as  a  redistribution  over  the  range  of  frequencies. 
The  ’shift’  in  the  acoustic  spectrum  is  further  highlighted  in  Figure  12,  demonstrating  the  contribution  of 
each  microphone  location  over  a  range  of  frequencies  from  Stp  =  0.2,  where  the  spectra  peak,  to  Sto= 3.0, 
where  the  magnitude  diminishes.  More  clearly  seen  here,  at  low  frequency  (Stp  =  0.2),  the  hot  jet  exhibits 
a  larger  noise  contribution  over  all  microphones,  peaking  at  the  downstream  microphone  (0=15°).  At  high 
frequency  {Sto  =  3.0),  the  opposite  is  seen  as  the  cold  jet  contributes  more  noise,  at  all  microphones,  peaking 
at  0=75°.  In  the  range  of  frequencies  between  these  two  states,  a  region  of  transition  is  seen  as  a  marked 
decrease  in  magnitude  is  observed  at  the  downstream  microphone,  with  increasing  frequency.  At  the  lowest 
frequency  presented,  the  noise  level  is  a  maximum  for  both  jets,  with  larger  magnitudes  exhibited  in  hot  jet 
at  all  angles.  The  contribution  then  diminishes  with  increasing  frequency,  more  drastically  in  the  hot  jet, 
ultimately  falling  below  the  values  of  the  cold  jet  at  all  microphones  at  Sto =3.0.  An  effect  which  is  shown 
to  begin  at  the  0=15°  microphone  location  and  cascade  to  consecutive  locations  (increasing  0  separation). 
The  AdD  difference  between  Sto= 0.2  and  0.3  is  approximately  35 dB  in  the  hot  jet,  at  0=15°,  and  near 
29 dB  in  the  cold  jet.  The  largest  aft  angle,  0=90°,  exhibiting  only  a  14 dB  and  11  dB  change  over  similar 
frequency  range,  in  hot  and  cold  jet,  respectively. 

A  striking  result  is  the  similar  distribution  over  all  microphone  locations  seen  for  both  the  cold  an  hot 
jet,  as  similar  trends  arc  illustrated  in  each  set  of  curves.  Where  suggested  previously  that  only  shallow 
aft  angles  experience  a  redistribution  of  acoustic  energy,  it  is  apparent  here  that  all  microphone  locations 
exhibit  an  increase  in  the  amount  of  low  frequency  noise  detected  and  a  reduction  in  high  frequency  noise,  in 
the  hot  jet.  Consistent  with  the  amplitude  changes  seen  at  low  frequency  in  the  spectral  densities,  a  larger 
increase  in  noise  level  at  low  frequency  is  demonstrated  at  shallow  aft  angles.  The  largest  low  frequency 
difference  between  the  cold  and  hot  jet  is  seen  at  Sto =0.2  and  0=30°,  of  dB.  Also  notable,  however 
not  as  clearly  disccrnable  from  the  spectra,  is  the  larger  decrease  in  high  frequency  noise  at  large  aft  angle, 
again  ~4 dB  seen  at  Sto= 3  0  at  0=75°.  This  coupled  with  the  more  pronounced  changes  in  amplitude  with 
increasing  frequency  seen  at  shallow  aft  angles,  substantiates  the  notion  that  these  angles  arc  dominated  by 
low  frequency  noise  contribution.  At  high  frequency,  the  higher  noise  levels  seen  at  large  aft  angles,  suggests 
that  the  opposite  is  also  an  accurate  assumption. 

An  examination  of  the  overall  sound  pressure  level  of  both  jets,  shown  in  Figure  13,  reveals  that  not  only 
is  the  acoustic  energy  redistributed  to  different  frequencies,  but  there  exists  an  overall  increase  in  the  level 
of  noise  produced  by  the  hot  jet  at  nearly  all  observed  angles,  when  integrated  over  frequency. 

SPL(f)  =  101og10  08) 

OASPL  =  101ogIO  (19) 

where, 

Prms  ~  Pvar  ~  J  H mm{f)4f  (20) 

This  trend  is  consistent  with  experimental  findings  of  Viswanathan27  and  Tanna,38  and  computational  work 
of  Bodony  k  Lclc.56  The  largest  increases  of  approximately  2 dB  is  seen  at  shallow  aft  angle,  the  0=30° 
microphone  location  exhibiting  a  maximum  AdB,  falling  nearly  linearly  with  increasing  separation  to  0=90°, 
where  a  slight  decrease  of  nearly  half  a  dB  is  observed.  Again  affirming  that  the  low  frequency  noise  is 
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(a) 


(b) 

Figure  11.  (a)Normalized  far-field  pressure  spectral  density  of  microphones  positioned  at  <f>=15°  &  90°  (b) 
and  sound  pressure  level  spectral  density. 
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Figure  12.  Sound  pressure  level  distribution  over  all  microphone  locations  as  a  function  of  Sto- 
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Figure  13.  Overall  sound  pressure  level  distribution  of  microphone  locations. 
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prevalent  at  shallow  aft  angles  and  high  frequency,  at  large  aft  angles,  as  the  high  frequency  noise  decrease 
in  the  hot  jet  is  depicted  here  in  the  0=90°  microphone  location. 

Tam  et  al  23  has  attributed  the  directional  distribution  of  the  axisymmctric  jet’s  far-ficld  acoustic  spec¬ 
trum  to  a  contribution  from  two  independent  components  of  the  flow-field,  of  varying  length-scales  and 
time-scales.  Although  jet’s  far-ficld  spectrum  is  thought  to  be  governed  by  the  contribution  of  large-scale 
structures,  Tam  &:  Auriault57  established  the  ’fine’  or  small  scales  (FS),  importance  to  the  weaker  and  more 
broadband  distribution  seen  at  large  aft  angles,  having  fairly  uniform  directivity.  The  larger  structures  (LS) 
arc  shown  to  contribute  to  the  dominant,  more  narrow- band  spectra  seen  at  shallow  aft  angles.  Universal 
similarity  noise  spectrum  for  each  component  empirically  derived  from  an  examination  over  a  wide  range  of 
databases,  has  demonstrated  similar  results  as  experimental  measurements,  Viswanathan;27  the  LS  spectrum 
fitting  measured  spectra  radiated  in  the  downstream  direction  at  shallow  aft  angle  and  the  FS  spectrum  fit¬ 
ting  that  observed  at  large  aft  angle.  As  a  hypothesis  for  this,  it  proposed  that  the  rapid  non-linear  decay  of 
instability  waves  with  the  collapse  of  the  potential  core  region,  result  in  an  intense  Mach  wave-like  radiation 
in  the  downstream  direction  (wavy- wall  theory).  As  the  dominant  structures  in  this  region  are  large  sealed, 
it  is  concluded  that  it  is  predominantly  their  contribution  which  is  responsible  for  the  observed  noise  in  this 
direction. 

However,  currently  there  is  no  widely  accepted  theory  to  interpret  the  mechanism  responsible  for  the 
change  in  spectral  shape  at  shallow  and  large  aft  angles.  Objections  coming  from  Ribner58  and  Fisher,59 
bolster  first  principles  of  jet  noise  theory,  emphasizing  that  this  pattern  is  based  on  refraction  effects  of  mean 
velocity  gradients.  The  predicted  ’figure-eight’  pattern  elongated  by  convective  amplification,  resulting  in 
the  dominant  distribution  seen  at  shallow  aft  angles.  Fisher59  points  out  the  lack  of  consistency  between 
the  universal  FS  and  LS  spectra  with  hot  jets  at  low  acoustic  Mach  number,  where  the  spectra  shows  an 
decrease  in  noise  levels,  rather  than  an  increase.  This  is  also  echoed  by  Viswanathan,41  who  although  not 
a  supporter  of  the  intrinsic  mechanism  driving  the  shifting  of  frequencies  in  hot  jets,  shows  discrepancy 
between  measured  data  and  LS  similarity  spectra  for  a  low  Mach  number  jet  over  a  range  of  temperature 
ratios. 

At  the  present  Ma=0.6  study,  there  is  however  evidence  seemingly  supporting  the  FS  and  LS  theory. 
The  spectra  in  both  cold  and  hot  jet  arc  of  similar  shape,  with  a  pronounced  increase  in  low  frequency 
and  decrease  in  high  frequency  contribution  observed  in  the  hot  jet.  The  properties  of  the  hot  jet  arc  such 
that  increased  mixing,  due  to  increased  spatial  growth  rate  of  instability  waves,  results  in  an  increased 
contribution  of  large  scaled  structures.  The  shift  to  a  dominant  low  frequency  column  mode  observed  in  the 
near-field  pressure  signature,  along  with  the  increase  in  low  frequency  acoustic  energy  seen  predominantly  at 
shallow  aft  angles,  may  serves  as  evidence  of  the  the  large  scale’s  contribution  predominantly  to  the  shallow 
aft  angles.  This  will  be  examined  further  in  the  next  section,  coupling  the  near  and  far-field  pressure,  as 
well  as  in  the  analysis  of  the  source  velocity  field. 


F.  Near-field  and  acoustic  far-field  coupling 

Although  the  near-field  azimuthal  array  does  not  lie  within  the  well  establish  ’noise  source  region’  of  6-10D 
(at  x/D= 2  and  r/D  ~  1),  the  pressure  signature  may  still  be  used  to  garner  information  about  the  acoustic 
far-field.  An  analysis  is  now  undertaken  to  examine  the  relationship  between  the  two  fields,  with  the  aim  of 
extracting  distinctions  in  dominant  frequency  of  the  structure  of  each  jet. 

The  two  signals  exhibit  weak  coherence,  at  a  maximum  of  less  than  10%  with  the  in-plane  transducer 
(0=90°)  and  the  downstream  microphone  6,  at  0=15°,  in  both  the  cold  and  hot  jet.  As  a  result,  coherency 
plots  arc  not  presented  here.  The  magnitude  of  the  cross-spectral  densities,  instead,  are  evaluated  using  only 
the  in-plane  transducer  and  microphones  1  &;  6  (0=90°  &  15°,  respectively).  These  are  calculated  as  follows, 
(k=l-6,  denotes  microphone  location,  j=l,  denote  in-plane  transducer) 

Gmp(/)  =  l;<P*(/)p*(/)>  (21) 

The  frequencies  over  which  these  signals  arc  coherent,  intuitively  should  coincide  with  those  exhibited  over 
the  range  of  k~ 2  slope  seen  in  both  the  near-field  spectra  (Fig.  6)  and  far-ficld  spectra  (Fig.  10),  indicating 
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Figure  14.  Near-field  and  acoustic  far-field  pressure  cross-spectra  ( single  in-plane  transducer).  (a)Cold  jet. 
(b)Hot  jet. 


an  acoustic  signature.  This  is  not  the  case,  as  the  amplitude  of  the  cross- spectra,  shown  in  Figure  14, 
demonstrates  a  pronounced  peak  over  a  range  of  Sto=  0.2-0. 5.  These  peaks  corresponding  instead  to  the 
range  of  frequencies  spanning  the  peak  energy  containing  regions  of  each  signal  (near-field  and  far-field 
spectra).  However  an  examination  of  the  extremely  low  magnitudes  confirms  trends  in  the  coherency  spectra 
alluded  to  earlier,  of  the  weak  correlation  between  the  two  signals.  A  10%  filter  is  applied  to  the  curves  for 
clearer  presentation.  It  is  interesting  to  note  that  the  correlated  energy,  however  limited  in  magnitude,  is 
a  result  of  a  predominantly  hydrodynamic  near-field  signature  with  the  acoustic  far-field.  The  existence  of 
this  correlation  suggests  that  the  pressure  measured  here  exhibits  some  signature  of  the  acoustic  source.  A 
similar  result  has  recently  been  shown  by  Hall  et  al 52 

In  comparing  the  two  plots  of  the  cold  and  hot  jet,  a  more  broadband  signature  is  again  seen  in  the  cold, 
exhibiting  a  weaker  correlation.  The  hot  jet  demonstrates  a  more  narrow  distribution  over  frequency,  falling 
to  negligible  magnitude  near  Sto=0A,  where  the  cold  extends  to  a  higher  value  of  Sto= 0.6  (as  compared  to 
peak  energy).  Of  particular  importance  is  the  increase  in  peak  frequency  with  increasing  aft  angle  separation, 
seen  in  both  jets,  in  concert  with  a  decrease  in  correlated  energy.  As  results  shown  in  the  analysis  of  the 
far-field  pressure  indicate,  the  dominant  frequency  of  the  acoustic  signature  also  increases  with  increasing 
aft  angle  separation.  The  peak  exhibited  here  at  <p= 15°  (mic  6),  coincides  with  the  dominant  low-frequency 
component  exhibited  in  the  near-field  pressure  signal,  shown  earlier.  This  trend  suggests  a  relationship 
with  the  energy  content  of  the  column  mode  of  the  near-field  pressure.  The  hot  jet,  exhibiting  a  greater 
correlation  here,  has  demonstrated  a  larger  amount  of  energy  in  mode  0  ( dominant  modal  contribution ), 
where  its  contribution  in  the  cold  jet  is  diminished  ( similar  modal  contribution  as  the  helical  mode),  and 
exhibits  lesser  correlation. 

To  investigate  this,  a  modal  decomposition  is  performed  in  similar  fashion  as  in  the  analysis  of  the 
near-field  pressure,  using  the  real-part  of  thccross-spectra,  Fourier  transformed  in  azimuth  (A0  — ►  m).  All 
transducer  locations  in  6 ,  and  a  single  microphone  location  in  0,  arc  used  to  realize  an  azimuthal  modal 
distribution,  as  a  function  of  each  microphone  location.  The  first  four  modes  arc  displayed  for  both  the  cold 
and  hot  jet,  again.  The  trend  previously  observed  in  the  cross-spectra  is  echoed  in  the  modal  distribution, 
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with  the  cold  jet  demonstrating  correlation  over  larger  frequency  range  and  the  hot  jet  showing  larger 
magnitudes  at  both  microphone  locations.  A  shift  to  lower  frequency  with  increasing  aft  angle  separation  is 
again  seen  in  both  jets,  as  the  peaks  at  microphone  1  (0=90°),  occur  at  higher  frequency,  near  Stp=0A,  and 
fall  to  Sto= 0.2  at  microphone  6  (0=15°).  The  correlated  energy  is  also  consistently  greater  at  microphone 
6.  In  examining  the  modal  contribution  across  the  two  studies,  the  mode  0  contribution  in  the  hot  jet 
at  microphone  6,  demonstrates  a  considerably  greater  increase  in  amplitude,  nearly  quadrupling  it’s  cold 
counterpart,  where  as  at  microphone  1  (0=90),  a  doubling  effect  is  seen.  This  is  consistent  with  the  increase 
in  low  frequency  noise  seen  at  all  microphone  locations  in  the  far-field  acoustic  spectra,  more  prominently 
seen  however,  at  shallow  aft  angle.  The  magnitude  of  all  others  modes  remains  nearly  consistent  for  both 
jets.  The  decrease  in  higher  frequency  energy  at  larger  aft  angles  in  the  hot  jet  acoustic  spectrum,  is  not 
clearly  seen  here,  though  it  can  be  inferred  from  the  trends  of  the  curves.  At  mic  6,  the  curves  extend  to 
higher  frequency  in  the  cold  jet  (Sto= 1.3)  than  in  the  hot  jet  (Sto =0.9).  A  more  pronounced  change  may 
be  evident,  were  it  not  for  the  lack  of  correlated  energy  at  the  higher  frequencies  exhibited  in  the  acoustic 
far-field  Stp  >1.  Recalling  that  the  decrease  in  high  frequency  noise  contribution  was  not  clearly  exhibited 
until  Sto= 3  in  Figure  12. 

The  directional  dependence  of  the  far-field  spectrum  with  frequency  is  highlighted  in  this  analysis.  It 
is  also  shown  that  it  is  the  low-frequency,  large-scale  column  mode  0,  prevalent  in  the  near-field  signature, 
that  drives  the  correlation  with  the  far-field,  especially  at  shallow  aft  angles.  (Trends  behave  similarly  to  the 
downstream  pressure  and  acoustic  field  coupling  seen  in  Hall  et  al.52).  At  a  maximum  correlation  with  the 
downstream  microphone  positioned  at  0=15°,  it  is  suggested  that  the  increased  contribution  of  mode  0  seen 
in  the  hot  jet,  is  responsible  for  the  increase  in  low  frequency  acoustic  energy  observed  at  shallow  aft  angles. 
At  increasing  aft  angles,  the  correlation  between  the  signals  decreases,  as  these  angles  have  been  linked  to 
the  high  frequency  contribution  of  the  ’fine-scale’  turbulence  structures.  Conversely,  these  results  imply  that 
it  is  the  higher  ordered  modal  structures  that  contribute  to  high  frequency  noise. 

The  traditional  quadroplc  distribution  of  jet  noise  docs  not  account  for  the  fluctuating  entropy  source 
term  (p  -  Oqp),  as  modeled  by  Morfey  et  al.40  to  exhibit  a  dipole-like  distribution.  Its  demonstrated 
contribution  of  low-frequency  energy  to  the  acoustic  spectrum  makes  this  a  viable  candidate  in  accounting 
for  the  shift  in  the  acoustic  spectrum  to  low  frequency,  and  the  overall  increase  in  noise  level  of  low  Mach 
number  heated  jets.  Where  the  reduction  of  the  overall  sound  pressure  level  seen  at  Ma  >0.7  is  concerned, 
much  research  is  left  to  be  done. 


G.  Source  velocity  field 

The  fluctuating  velocity  fields  of  both  the  cold  and  hot  jet  arc  now  compared  in  an  analysis  aimed  at  iden¬ 
tifying  the  dominant  features,  intrinsic  to  each,  that  contribute  to  sources  of  noise  .  The  use  of  stereoscopic 
PIV  to  sample  the  velocity  field  allows  superior  spatial  resolution  as  compared  to  other  tools.  However 
it’s  low  sampling  rate,  makes  acquisition  of  temporal  time  scales,  in  high  Reynolds  number  flows,  difficult. 
Thus  the  sampling  rate  used  (1  Hz),  was  chosen  to  ensure  statistically  independent  snapshots,  rather  than 
attempting  to  resolve  dominant  time  scales. 

1.  Single-point  statistics 

The  mean  velocity  profiles  at  all  measured  locations  (x/D=4.5-10)  are  shown  in  Figure  16a,  comprising  an 
average  over  a  minimum  of  1980  vector  snapshots.  The  average  centerline  velocity  at  the  furthest  upstream 
location  is  reported  at  204m/s  for  the  cold  jet  and  199m/s  in  the  hot  jet.  As  the  present  study  fixes  acoustic 
Mach  number,  and  thus  exit  velocity,  these  values  are  in  concert  with  expected  the  Uj  =  A/a/aoo=204m/s 
at  A/a=0.6.  Although  the  average  velocity  exhibited  in  the  hot  jet  is  well  within  the  uncertainty  of  the 
measurement,  it  is  now  assumed  that  this  value  corresponds  to  a  premature  centerline  decay,  a  known 
feature  of  hot  jets.  Due  to  limitations  in  optical  access  near  the  jet  exit,  experimental  verification  was  not 
possible.  However,  previous  measurement  using  LDA  and  PIV  have  reported  centerline  values  of  203 m/s  and 
202m/s,  respectively.  The  following  presentation  of  the  mean  statistics  follows  in  this  assumption,  whereby 
the  centerline  velocity  of  the  cold  jet  (U^°ld^),  is  utilized  for  all  normalization. 
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(b) 

Figure  15.  Modal  distribution  of  near-field  pressure  and  acoustic  far-field  cross-spectra.  (a)Cold  jet.  (b)Hot 
jet. 
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(b) 

Figure  16.  (a)Mean  velocity  distribution  of  (left)hot  jet  and  (right)cold  jet.  (b)Mean  velocity  distribution  of 
both  jets  at  (left)x/D=4.5  and  (right)x/D=10. 
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Apparent  in  the  profiles  is  the  more  rapid  decrease  in  centerline  velocity  with  downstream  distance,  seen 
in  the  hot  jet,  indicative  of  an  increased  decay  rate,  as  seen  by  Tanna,38  Bridges  &  Wcrnct,60  and  Lew  et  a/.61 
Recently,  Panda  et  ai72  showed  through  measurement,  that  with  the  addition  of  heat,  the  axisymmetric  jet 
exhibits  increased  density  fluctuations,  which  is  linked  to  an  increased  growth/spreading  rate.  This  is  not 
seen  here,  as  both  sets  of  curves  extend  to  similar  radial  locations  at  similar  strcamwisc  location.  From  the 
curves,  a  spreading  rate  near  7.5°,  (0.13x)  is  calculated.  However,  the  rapid  decay  of  the  hot  jet  is  shown 
in  Figure  16b  to  decrease  the  half- width  of  the  curves  near  their  apex,  and  as  such  produces  a  larger  inner 
shear  layer  ( greater  potential  core  collapse  angle).  Ultimately,  a  larger  mixing  layer  is  exhibited  in  the  hot 
jet. 

The  trend  observed  indicating  a  shorter  potential  core  length,  is  investigated  further  as  the  centerline 
velocity  of  each  jet  is  presented  as  a  function  of  streamwise  location  in  Figure  17a  and  tabulated  in  the 
Table  4  below,  recalling  again  that  Uci  refers  to  that  exhibited  by  the  cold  jet  at  x/D=4-5.  With  the  criteria 


Table  4.  Centerline  Velocity 


x/D 

%Uci  Cold  Jet 

%Uci  Hot  Jet 

4.5 

100 

97.41 

5 

97.94 

94.49 

5.5 

97.20 

92.79 

6 

94.12 

88.76 

of  0.95 Uci  indicating  the  collapse  of  the  potential  core  region,  the  curves  indicate  that  the  hot  jet  collapses 
further  upstream,  near  x/D= 5,  while  the  cold  jet  is  seen  to  persist  until  x/D= 6. 


The  Witze62  correlation  parameter  is  an  effective  means  of  comparing  jets  of  varying  temperature  and 
Mach  number,  minimizing  compressibility  effects,  defined  by, 


W  =  2/c 


/  X  -  Xc\  /  £oo  \  ^ 

l  D  J\Pi) 


(22) 


where 


k.  =  0.08(1  -0.16*  Ud/a^y0-22  (23) 

Here  xc  refers  to  potential  core  length  defined  as  axial  location  of  0.95f/c/,  pj  and  poo  are  jet  and  ambient 
density  respectively,  doo  the  ambient  sound  speed,  and  D  is  jet  diameter.  The  curves  collapse  rather  well, 
not  surprisingly,  considering  the  fairly  low  temperature  ratio  and  small  range  over  which  measurements  arc 
acquired.  The  streamwise  velocity  curves  are  also  plotted  using  the  similarity  variable  t](x)  =  (r  -  ro.s)/x, 
where  r0 .5  marks  the  location  of  0.5 Uci,  in  Figure  17b.  The  profiles  for  each  jet  collapses  well,  onto  a  similar 
curve.  The  region  of  disparity  seen  between  ?7(x)=-0.1  to  -0.05  indicates  a  lost  of  similarity  consistent  with 
the  collapsing  of  the  potential  core  region.  The  collapse  of  the  curves  using  these  various  scaling,  shows  the 
data  to  be  accurate  to  first  order. 

The  turbulence  intensities  for  both  the  cold  and  hot  jet  are  presented  in  Figures  18a  &  b,  as  a  function 
of  radial  distance  r/D  and  rj(x).  Streamwise  locations  of  x/D= 4.5,  6  &  8,  are  chosen,  corresponding  to  the 
existence,  collapse,  and  fully  mixed  potential  core  region.  The  top  curves  corresponding  to  the  streamwise 
component  (u),  middle,  the  radial  (v),  and  bottom,  tangential/azimuthal  (w).  When  referenced  against 
previous  measurements,  these  values  fall  slightly  below  those  reported,  at  just  over  11%  for  the  streamwise 
component  and  6%  for  the  in-plane  velocities  (16%  streamwise,  7%  radial  and  azimuthal  from  LDA,  and  15% 
streamwise,  7%  radial  and  azimuthal  from  PIV).  There  is  an  increase  in  turbulence  levels  seen  in  the  hot  jet 
as  the  streamwise  component  increases  to  nearly  13%,  with  a  slight  increase  in  magnitude  seen  in  the  other 
velocity  components  38-60-61*.56  The  increased  mixing  of  the  hot  jet  is  again  apparent  as  the  turbulence 
levels  at  the  centerline  no  longer  fall  at  xjD— 8,  as  in  the  cold  jet,  signifying  a  completely  mixed  region. 
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(b) 


Figure  17.  (a)Mean  centerline  velocity  decay  (left)x/D  scaling  and  (right)Witze  scaling.  (b)Collapsed  stream- 
wise  velocity  profiles  of  (left)cold  jet  and  (right)hot  jet,  using  ro.5  scaling. 
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The  diminished  magnitudes  arc  coupled  with  evidence  of  pronounced  turbulent  fluctuations  along  the 
edge  of  the  outer  shear  layer.  The  disturbance  is  confined  to  a  radial  location  along  the  low-speed  side  of 
the  mixing  layer  and  appears  to  grows  with  the  jet’s  evolution.  This  is  more  prominent  in  the  streamwise 
component,  indicating  the  presence  of  a  ’streamwise  disturbance.  It’s  influence  is  less  prominent  in  the  hot 
jet  where  an  increased  turbulence  intensity  is  exhibited.  Recalling  the  placement  of  the  transducer  array 
near  the  jet  exit  (x/D= 2,  r/D  c^l)  in  Figure  3,  it  is  clear  that  it’s  positioning  is  the  most  likely  culprit  of 
this  disturbance,  generating  streamwise  vortices  acting  similarly  to  a  chcvron/microjet  nozzle.  The  holders 
for  the  individual  transducer  probes  consist  of  small  aluminum  blocks  that  arc  held  in  place  by  two  thin,  flat 
sheet  metal  strips,  on  either  side.  These  are  screwed  together,  clamping  the  blocks  in  between.  These  thin 
strips  act  as  knife  edges,  perturbing  the  entrained  air  near  the  jet  lip. 

Arakeri  et  a/.64  have  reported  decreases  in  turbulence  intensity  with  the  use  of  microjets,  with  injection 
as  minimal  as  1%  of  mean  core  flow.  Alkislar  et  al .63  concluded  that  this  influence,  although  generating 
strong  streamwise  vorticity  initially,  decayed  rapidly  with  downstream  distance  confining,  it  to  a  smaller 
region  of  the  jet. In  the  present  ease,  the  array’s  positioning  ensures  there  is  no  penetration  into  the  core  flow 
of  the  jet,  resulting  only  in  a  perturbation  of  the  outer  edge  of  the  forming  shear  layer,  with  disturbances 
confined  to  this  region.  Therefore,  no  change  is  exhibited  in  the  mean  statistics  and  the  increased  growth 
rate  with  respect  to  the  natural  jet  is  also  not  seen  here,  as  comparisons  of  mean  velocity  profiles  and  shear 
layer  thickness  to  experiments  in  the  same  facility  at  similar  Mach  number  arc  nearly  completely  matched. 


2.  Reynolds  stresses 

With  the  predisposition  of  investigating  acoustic  source  features,  the  normalized  Reynolds  stresses  for  the 
axial  location  x/ D  =  8,  lying  within  the  noise  producing  region,  arc  presented  in  Figures  19  &:  20.  These 
are  computed  using, 


&j(r> r'.0) 


(Ui(r, 9,  t„)uj(r',  9  +  A0,  ts)) 

vl 


(24) 


where  ts  denotes  snapshot  time  step,  Uci,  the  jet  exit  centerline  velocity,  and  index  notation  is  adopted,  as 
(ui,U2,  and  U3)  replaces  the  previous  (u,v,w)  notation  ( corresponding  to  streamwise ,  radial,  and  tangential 
directions,  respectively).  The  <  .  >  signifies  an  average  over  snapshot  timing,  as  well  as  similar  A 9  separation. 


Although  the  full  Reynolds  stress  matrix  is  available,  only  the  radial  grid  representing  the  normal  stresses 
(£,,)  at  a  separation  of  A0=O,  are  presented  in  the  contours  in  Figure  19.  The  streamwise,  or  axial,  stresses 
arc  shown  to  dominate  the  flow  field,  with  the  radial  stresses  consistently  exhibiting  the  largest  spatial  range. 
The  peaks  in  the  curves  of  the  cold  jet  arc  centered  along  the  turbulent  mixing  layer  (r  =  r'  =  0.5),  while 
the  hot  jet  demonstrates  a  trend  consistent  with  a  more  rapid  decay  of  potential  core  region,  as  peaks  have 
moved  closer  toward  the  centerline  at  this  downstream  location.  Consistent  with  this,  the  length  scales  are 
shown  be  larger  in  the  hot  jet,  as  compared  to  it’s  cold  counterpart. 

At  an  azimuthal  separation  A0=O°,  along  the  line  r=r'  ( i.e .  the  diagonal  of  the  contours ),  similar  curves 
as  depicted  by  the  turbulence  intensities  presented  previously  can  be  realized.  These  are  presented  in  the 
set  of  curves  beneath  the  stress  contours.  (Note  that  magnitudes  of  the  y-axis,  although  similar,  are  of 
different  scaling  for  a /Uci  and  £).  Paralleling  the  increased  turbulence  levels,  the  axial  stresses  in  the  hot 
jet  demonstrate  an  increased  contribution  with  the  addition  of  heat.  It  is  important  to  again  highlight  the 
similar  spreading  rates,  as  the  curves  reach  a  minimum  at  similar  values  along  the  horizontal  axis  in  both 
sets  of  curves.  The  shear  stress  components  are  presented  in  Figure  20,  illustrating  that  only  the  £12  term 
has  a  significant  amount  of  energy.  Similar  trends  arc  observed  here,  as  the  magnitude  of  the  stresses  and 
spatial  length  scales  are  again  greater  in  the  hot  jet. 

Where  noise  production  is  concerned,  the  Reynolds  stresses  can  be  viewed  as  a  simplified  representation 
of  the  Lighthill’s  stress  tensor,  =  pUiUj.  Although  not  a  valid  approximation  to  the  complete  source 
term  with  the  addition  of  heat,  the  weighting  of  these  terms  can  be  examined  to  give  a  measure  of  there 
contribution  to  overall  noise  production.  The  increased  magnitudes  consistently  seen  in  the  hot  jet  would 
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Figure  18.  Turbulence  intensity  profiles  (a)Cold  jet  (b)Hot  jet. 
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(b) 


Figure  19.  Normal  Reynolds  stresses  at  x/D  =  8  (a)Cold  jet  (b)Hot  jet. 
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r/D 


(b) 

Figure  20.  Shear  Reynolds  stresses  at  x/D  —  8  (a)Cold  jet  (b)Hot  jet. 
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seem  to  indicate  a  clear  link  between  these  stresses  and  the  increase  in  overall  sound  pressure  level  shown 
earlier  in  an  analysis  of  the  far-ficld  acoustics.  However,  as  this  term  scales  with  density,  the  vast  difference 
in  mean  density  in  each  jet  (pcoid= 1  428  and  phot  =0.759)  demonstrates  that  this  is  not  the  case.  As  such, 
the  Reynolds  stress  terms  contribution  to  noise  generation  in  the  hot  jet  would  be  diminished  as  compared 
to  the  cold  jet.  This  again  suggests  the  presence  of  an  additional  acoustic  source  with  the  addition  of  heat. 


3.  POD  analysis 

Known  in  mathematical  circles  as  Karhunen-Loeve  decomposition  or  principle  component  analysis ,  POD 
serves  a s  a  method  for  extracting  a  basis  set  of  modal  functions  from  an  ensemble  of  signals.  Where  the 
turbulence  community  is  concerned,  its  serves  to  bring  the  orderly  motion,  of  a  seemingly  random  flow 
field,  to  the  forefront.  Similar  to  an  energy  budget,  the  POD  greatly  reduces  the  dimension  of  turbulence 
by  decomposing  the  structure  as  a  function  of  overall  energy  contribution.  Based  on  the  spectral  theory 
of  compact,  adjoint  operators,  the  method  generates  a  basis  set  of  spatial  eigenfunctions  resulting  from 
the  solution  to  the  Fredholm  eigenvalue  problem.  The  ensemble  of  signals  dictates  these  to  be  the  most 
appropriate  functions,  resulting  from  a  maximization  of  the  signals’  largest  mean  projection  onto  this  basis 
set. 


We  start  by  seeking  a  deterministic  set  of  functions  that  best  describes  our  ensemble  of  signals  (a  velocity 
field,  u(x).  From  this  set,  a  ’candidate  structure’  0(x),  which  maximizes  the  normalized  inner  product  with 
the  velocity  field,  is  found.  Vector  notation  is  dropped  for  simplicity  in  the  formulations  below. 


<](«■  •■/pi2)  _  (i(tt-<ft)i2) 
(0-0*)  (0'0*) 


(25) 


This  formulation  leads  to  the  Fredholm  integral  eigenvalue  problem,  as  a  necessary  condition  for  this  to  hold 
is  that,  0(x),  be  eigenfunctions  of  the  two-point  spatial  correlation  tensor  of  u(x). 


/  (ui(x,t)uj(x' ,t))  <pj(x)dx  =  \<pi(x)  (26) 

JD 

the  two-point  spatial  correlation  tensor  known  as  the  kernel  is  generally  replaced  by  Rlj(x^xf).  Hilbcrt- 
Schmidt  theory  dictates  that  as  a  result  of  the  inner  product  formulation  in  Hilbert  (complex  vector)  space, 
the  kernel  Rij  will  be  a  Hermitian  operator.  Additionally,  solutions  to  the  problem  will  consist  of  a  infinite 
set  of  eigenfunctions,  0(n)(x),  and  a  corresponding  ordered  sequence  of  eigenvalues  X^n\  where  n  denotes 
order  in  the  sequence.  The  properties  of  the  solution  arc  such  that, 

•  The  maximization  in  Eq.(  25)  corresponds  to  the  first  eigenvalue,  A^ 

•  The  set  of  eigenfunctions,  </>[n)  are  orthonormal, 


[  <t>^{x)^\x)d{x)  =  6pq 
JD 


(27) 


•  The  projection  of  the  original  ensemble  onto  the  set  of  eigenfunctions  yields  a  set  of  coefficients  an , 


an(t)  =  [  ut(f,t)0(n)(x)<ff 
JD 


(28) 


•  These  POD  coefficients  an  will  be  random  and  uncorrelatcd, 

(<Z  nflm)  =  A*  ^  &nm 


(29) 


•  The  original  ensemble  can  be  reconstructed  by  projecting  the  set  of  random  POD  coefficient  an  onto 
the  set  of  eigenfunctions, 


Ui(x,t)  =  ^an(t)0Sn)(x)  (30) 

n=  1 
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The  dimension  of  the  kernel  TV,  in  the  eigenvalue  problem  is  equivalent  to  the  size  of  the  domain  (£), 
which  can  be  considerable  large.  The  present  formulation  of  a  velocity  field  u(x)  would  triple  in  dimension 
if  each  component  was  included,  (nc  •  TV).  Sirovich31  realizing  the  inherent  difficulties  from  a  practical  point 
of  view  attempted  to  derive  a  less  computationally  costly  solution,  which  he  referred  to  as  ’snapshot’  POD. 
The  general  basis  of  the  problem  is  the  same  as  in  the  ’classical’  form,  with  the  apparent  distinction  being 
the  formulation  of  the  kernel,  and  subsequent  domain  over  which  the  eigenvalue  problem  is  solved.  We  again, 
begin  with  an  ensemble  of  a  velocity  field,  however  now  separated  by  some  time  scale  r,  on  the  order  of  or 
greater  than  the  correlation  time  such  that, 

it(m)(£)  =  u(£,  tut)  (31) 


Each  m  referring  to  a  ’snapshot’  or  sufficient  temporal  separation  making  the  velocity  field  uncorrclated 
in  m.  In  similar  fashion,  the  maximization  leads  to  the  integral  eigenvalue  problem,  with  the  kernel  now 
replaced  by  two-point  temporal  correlation  tensor  of  i^m)(£). 


1  M 

C(x,f')  =  -  ^  u('n)(x)u(m)(f') 

m—  1 

equivalently  written  as, 

C(t,t')  =  Tf.  I  Ui{x,t)Ui(x,t')dx 

1  Jd 


(32) 


(33) 


where  T  is  the  period,  D  is  the  domain  dictated  by  (z),  and  t'  represents  a  temporal  separation.  The 
resulting  integral  eigenvalue  problem  now  becomes, 


L 


C(t,  t')an(t')dt'  =  A<n)an(t) 


(34) 


The  eigenvalue  problem  is  now  solved  over  the  temporal  domain,  for  the  random  POD  coefficients,  an(t), 
now  more  clearly  described  as  temporal  eigenfunctions.  Some  key  aspects  of  the  solution  arc  highlighted 
below. 


•  The  POD  coefficients  an  will  again  be  random  and  uncorrelated,  as  a  property  of  the  basis  of  formu¬ 
lation. 


<anam)r  =  A(n)(5 


nm 


(35) 


•  The  projection  of  the  original  ensemble  onto  the  set  of  temporal  eigenfunctions,  an,  now  yields  a  set 
of  spatial  eigenfunctions  <t>[n\ 

4>|n)(z)  =  ifXbO  jTan{t)ui{x,t)dt  (36) 

•  The  coefficient  must  be  applied  to  satisfy  the  condition  that  the  set  of  spatial  eigenfunctions, 

(j)[n^  are  orthonormal, 

[  <t>\p){x)<p\v){x)d(x)  =  <5P,  (37) 

Jd 


The  snapshot  method  aims  to  reduce  the  dimension  of  the  integral  eigenvalue  problem  from  nc  ■  TV  to  T. 
However,  the  pronounced  difference  shown  in  the  formulation  of  the  two  techniques  may  ultimately  serve  to 
give  no  clear  advantage.  Both  techniques  can  be  applied  to  similar  data  sets.  Should  the  flow  configuration 
involve  a  coarser  spatial  grid  size  than  temporal  (TV  «:  T),  or  a  temporal  grid  that  is  significantly  dense, 
however  still  less  coarse  than  the  spatial  (T  <  TV,  TV  >  1),  then  the  snaphot  method  loses  its  advantage. 

The  snapshot  proper  orthogonal  decomposition  is  performed  giving  a  modal  decomposition  of  the  flow 
field,  identifying  the  most  energetic  or  ’coherent’  structures,  in  an  attempt  to  ascertain  a  measure  of  the 
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Mode  Number 


Figure  21.  (a)Convergence  of  POD  modes  (left)  at  x/D= 4.5  (right)  and  x/D= 8. 


large  scale  motion  of  the  flow- field;  as  it  is  the  contribution  of  these  structures  that  appears  to  be  the  marked 
difference  in  the  cold  and  hot  jet  The  ’snapshot’  version  of  the  method  is  used,  comprising  an  average  of 
1980  PIV  vector  maps  in  each  instance. 

The  formulation  of  the  kernel  matrix  in  the  POD  eigenvalue  problem  using  the  two-point  fluctuating 
velocity  correlations,  dictates  that  the  summation  of  the  total  number  of  eigenvalues  generated  gives  the 
total  turbulent  kinetic  energy  of  the  flow  field.  The  distribution  of  the  system’s  energy  in  each  jet,  is 
very  telling  as  to  the  dimensional  nature.  The  ordering  of  the  eigenvalues  requires  that  the  most  energetic 
structures  of  the  flow  field  will  be  represented  by  the  first  POD  mode.  In  Figure  21,  the  distribution  of 
energy  is  presented  for  both  jets  at  streamwise  locations  x/D= 4.5  &:  8.  The  hot  jet  is  shown  to  converge 
more  rapidly  at  both  locations,  depicting  a  higher  amount  of  the  system’s  total  energy  over  the  lower  modes. 
The  largest  disparity  is  seen  at  x/D=4.5,  where  the  first  mode  of  the  hot  jet  nearly  doubles  the  contribution 
of  the  cold,  exhibiting  more  than  13%  of  the  flow  field’s  total  energy.  A  similar  trend  is  seen  downstream, 
with  a  less  pronounced  difference.  The  turbulent  kinetic  energy  distribution  of  each  jet  is  tabulated  below. 
Percentages  represent  %  of  total  energy,  and  values,  the  number  of  modes.  The  low-dimensional  nature, 
characteristic  of  the  axisymmetric  jet  is  seen  here  as  90%  of  the  flow  field’s  total  energy  can  be  recovered 
using  an  average  of  less  than  7%  of  the  total  modes  in  all  jets  (1980  total  as  dictated  by  the  snapshot  method). 
Where  the  hot  jet  is  concerned,  it  is  seen  that  the  features  of  this  jet  are  ’more’  low-dimensional  in  nature 
than  that  of  the  cold, exhibiting  a  steeper  convergence  slope.  This  is  again  more  evident  at  the  streamwise 
location  x/D= 4.5,  closer  to  the  jet  exit.  In  the  next  section,  the  structure  of  these  ’modes’  is  evaluated  to 
gain  some  insight  into  the  ’physical’  flow  features  responsible  for  this  distribution.  Previous  references  to 
the  modal  structure  of  the  source  velocity  field  have  been  in  the  context  of  the  Fourier- azimuthal  modes.  In 
the  next  section,  the  flow  field  is  examined  as  a  function  of  the  POD  modal  structure.  The  ’structure’  of  the 
individual  POD  spatial  modes,  4 ^  refers  to  an  eigenfunction  representation  of  the  energy  in  the  flow  field, 
and  ultimately  not  the  jet’s  ’flow  structure’.  Interpretation  of  these  is  at  times  ambiguous.  The  features 
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Table  5.  %  Total  Energy  Convergence  ( Number  of  POD  modes) 


Jet/Location 

30% 

50% 

70% 

90% 

Cold,  x/D= 4.5 

7 

17 

37 

101 

Hot,  x/D= 4.5 

5 

12 

30 

93 

Cold,  x/D= 8 

8 

20 

48 

139 

Hot,  x/D= 8 

7 

17 

42 

131 

of  the  flow  will  also  be  evaluated,  by  summing  over  the  projection  of  the  POD  temporal  modes,  an(t)  onto 
these  eigenfunctions,  to  reconstruct  the  original  flow  field. 

oo 

Ui(x,  t)  =  ^2  an{t)4>[n) (x)  (38) 

n=l 

The  total  number  of  modes  used  in  this  projection,  dictates  the  total  amount  of  energy  in  the  representation 
of  the  flow  field.  As  such,  this  technique  allows  the  realization  of  the  most  energetic  ’flow  structure’  within 
the  random  flow  field,  using  only  the  1st  POD  mode,  greatly  reducing  the  order  of  the  system. 


4.  Modal  structure 

It  has  been  suggested  that  the  use  of  the  snapshot  POD  method,  along  with  reducing  computational  dimen¬ 
sions,  does  not  ’presume’  periodicity  of  the  flow,  as  associated  with  a  Fourier  decomposition.  It  follows  that 
the  ’snapshot’  POD  decomposition  gives  a  ’true’  representation  of  the  measured  jet’s  characteristics,  in  the 
sense  that  it  does  not  force  fourier  modes,  and  any  asymmetries  ( from  flow  field  or  measurement ),  would 
be  highlighted.  In  the  azimuthal  direction  ( 0 ),  if  the  measured  flow  field  is  truly  periodic  and  stationary 
in  time,  by  definition,  the  POD  decomposition  mathematically  reduces  to  a  Fourier  decomposition  in  this 
direction.  Should  the  symmetry  in  measurement  hold,  then  similar  structures  can  be  garnered  in  both  cases. 
However,  the  ’snapshot’  POD  method  used  couples  the  decomposition  of  both  the  radial  (r)  and  azimuthal 
direction  ( 9 ).  As  such,  care  must  be  taken  where  direct  comparisons  between  the  two  modal  structures  is 
concerned. 

Figures  22-  25  illustrate  the  flow  features  of  the  first  9  POD  modes  of  the  axial  velocity  component, 
<t>[n\  for  both  the  cold  and  hot  jet,  again  at  x/D= 4.5  and  8.  The  cold  jet  exhibits  a  dominant  structure 
resembling  a  Fourier-azimuthal  mode  1  structure  (  22a).  This  is  the  most  likely  structure  present  in  the  flow, 
as  it  represents  POD  mode  1,  contributing  the  most  energy  to  the  flow  field.  POD  modes  7  &  8  are  similar 
to  a  Fourier  mode  4  structure,  and  suggesting  a  phasing  relationship  of  approximately  45°  between  the  two. 
The  higher  mode  9  structure  is  similar  to  higher  Fourier-azimuthal  mode  5,  suggesting  an  ordering  similar 
to  that  of  the  Fourier  modes.  In  line  with  this  reasoning,  combinations  of  POD  modes  2-6  are  examined  to 
determine  if  lower  ordered  Fourier  like  modes  may  be  present  within  these  structures.  Modes  4  &  5  are  shown 
to  form  a  POD  mode  7-like  structure  (Figure  22b).  The  recent  analysis  of  Pinier67  highlights  this  result. 
Not  clear  here  is  whether  these  structures  result  from  the  radial  inhomogeneities  or  the  azimuthal.  The 
higher  ordered  Fourier-modal  distribution  over  POD  modes  present  in  1,13,33  is  not  seen  here  as  azimuthal 
directional  decomposition  is  not  performed  prior  to  application  of  the  POD. 
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Figure  22.  (a)POD  eigenfunctions  of  the  axial  velocity  component  <t>[n ^  ( 1st  9  modes),  of  the  cold  jet  at  xf D= 4.5. 
(b)Sum  of  4-  <t>[5^ . 
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Figure  25.  POD  eigenfunctions  of  the  axial  velocity  component  <t>^  ( 1st  9  modea),  of  the  cold  jet  at  x/D=  8. 

The  more  striking  result  is  the  clearly  well-organized  structures  present  in  the  hot  jet  at  this  streamwise 
location,  closer  to  the  jet  exit,  in  Figure  23.  The  only  flow  structures  exhibited  here  are  those  similar  to 
Fourier-azimuthal  mode  structures.  POD  mode  1  resembles  a  Fourier  mode  0,  the  column  mode,  as  the  most 
likely  structure  in  the  hot  jet.  The  higher  modes  are  paired,  representing  similar  Fourier  modes,  each  with 
a  characteristic  phasing  relationship.  Modes  2  &;  3  at  90°,  and  consecutively  higher  modes  at  a  decreasing 
phasing  angle.  This  ’anomaly’  is  hypothesized  next  as  the  streamwise  separation  is  increased  to  x/D= 8. 

At  x/D= 8,  the  structures  of  the  hot  jet  remain  fairly  well-organized,  exhibiting  a  larger  radial  extent 
(Figure  24).  The  dominant  Fourier-azimuthal  column  mode  has  now  switched  to  a  helical  mode  1,  also  seen 
at  POD  mode  2.  The  higher  modes  exhibit  similar  pairing  and  phasing,  as  seen  at  x/D= 4.5.  Mode  6  &  7 
however,  now  combine  to  form  an  azimuthal  mode  3  structure,  a  characteristic  seen  earlier  in  the  cold  jet. 
In  Figure  25  the  cold  jet  at  this  location  illustrates  a  behavior  similar  to  that  at  the  upstream  location,  with 
an  increased  radial  extent  and  the  presence  of  what  appears  to  be  higher  ordered  Fourier  modes  at  POD 
mode  9. 

The  observed  trend  now  allows  inferences  to  be  made  regarding  the  structure  of  the  hot  jet  previously 
shown  in  Figure  23a.  With  the  addition  of  heat,  at  a  fixed  acoustic  Mach  number,  a  lower  Reynolds  number 
is  effected.  As  the  flow  field  moves  to  a  more  ’laminar’  flow,  the  structures  become  more  organized  radially, 
resulting  in  a  more  homogeneous  and  periodic  jet,  as  the  axisymmetric  column  mode  structure  is  shown  to 
persist  longer  in  the  source  velocity  field  of  the  hot  jet.  This  is  in  line  with  the  previous  analysis  of  fluctuating 
pressure  field,  as  the  hot  jet  also  showed  a  pronounced  column  mode  contribution  as  compared  to  the  cold 
jet.  The  breakdown  of  the  organization  leads  to  the  trend  of  ’splitting’  Fourier  modes  observed,  where  the 
POD  modes  now  sum  to  form  Fourier  modes,  or  simply  exhibit  non-representable  eigenfunctions.  This  only 
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occurs  downstream  of  the  collapse  of  the  potential  core  in  the  hot  jet,  where  the  axisymmctric  shear  layers 
have  merged. 

5.  Velocity  field  reconstruction 

The  gain  a  better  insight  as  to  the  features  of  the  original  flow  field  that  arc  tied  to  the  representation  of 
the  POD  modal  structure,  the  axial  velocity  component  is  now  reconstructed  (Figures  26  &  27).  The  flow 
field  is  first  generated  using  a  single  mode  ( POD  mode  1,  n=l )  to  illustrate  more  clearly,  the  contribution  of 
the  most  energetic,  'most  likely ’,  flow  feature.  The  streamwise  location  of  x/D= 4.5  is  chosen,  as  this  region 
exhibits  the  largest  differences  between  both  jets. 
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Figure  26.  Velocity  field  reconstruction  of  axial  component  at  x/D= 4.5.,  in  the  cold  jet. 

The  colormap  scaling  is  normalized  to  the  range  of  minimum  to  maximum  magnitude  of  fluctuating  velocity 
at  this  snapshot  (±u’),  to  better  illustrate  the  overall  contribution  of  each  mode.  As  seen  in  the  figures, 
the  hot  jet  depicts  a  brighter  contour  indicating  a  larger  contribution  of  energy  in  the  first  POD  mode  (as 
dictated  by  the  steeper  convergence  curves  presented  in  Figure  21). 

This  trend  continues  until  the  50  mode  representation,  where  the  levels  of  the  contours  in  the  cold  jet 
begin  to  more  closely  resemble  that  of  the  original  velocity  snapshot.  This  level  of  similarity  is  reached  in 
the  hot  jet  at  n=5.  The  more  interesting  feature  here  is  the  greater  change  in  the  ’richness’  of  the  flow 
features  seen  in  the  cold  jet,  at  similar  mode  count.  The  original  velocity  field  in  both  jets  exhibit  a  complex 
structure,  with  peaks  lying  with  the  region  corresponding  to  the  mixing  layer.  In  the  center  of  the  curves,  the 
potential  core  is  represented.  The  hot  jet  depicts  a  smaller  core  region,  and  the  presence  of  larger  structure 
within  the  mixing  layer.  As  POD  mode  contribution  is  increased,  there  is  a  greater  disparity  between  the 
shapes  and  magnitudes  of  the  peaks  in  the  cold  jet  at  similar  mode  count,  as  compared  to  the  hot  jet.  This 
especially  apparent  at  lower  mode  count. 
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Figure  27.  Velocity  field  reconstruction  of  axial  component  at  x/D= 4.5.,  in  the  hot  jet. 


The  structure  of  the  hot  jet  remains  fairly  similar  over  the  addition  of  the  first  20  POD  modes,  depicting 
only  a  gradual  increase  in  the  contribution  of  total  energy.  The  energy  distribution  is  such  that  the  bulk 
flow  features  are  realized  using  a  lesser  number  of  modes.  The  implies  similar  scaling  of  the  low  modes, 
whereby  once  the  large  scale,  high  energy  contribution  is  realized,  the  higher  ordered  modes  contribute  little 
to  the  shape  of  the  curves,  until  near  n=20.  This  trend,  as  does  the  contours  of  the  original  velocity  map 
( although  only  an  instantaneous  snapshot ),  suggests  that  a  greater  disparity  in  the  scales  of  the  dominant 
flow  structure  exists  in  the  hot  jet  i.e.  mainly  large  scales,  as  the  low  energy  contribution  of  the  small  scales 
has  little  impact  on  the  flow  features.  The  cold  jet  conversely  exhibits  a  trend  of  a  greater  ’disparity’  in 
the  scales,  i.e.  ranging  in  scales  from  large  to  small  more  evenly  distributed,  exhibiting  more  pronounced 
changes  in  the  flow  features  with  the  addition  of  each  mode. 

Thus  far,  it  has  been  found  that  the  major  differences  between  the  hot  and  cold  jet  flow  structure  resides 
in  the  larger  length  scales  exhibited  in  the  hot  jet.  This  feature  is  again  highlighted  here.  The  modal 
structure  analysis  has  now  demonstrated  that  it  is  the  differences  in  the  dominant  modal  ordering  (mode  0, 
hot  and  mode  1,  cold)  of  the  two  jets  that  stands  as  the  source  of  this  scaling.  Larger  scales  persist  owing  to 
the  stagnation  of  the  natural  development  of  the  jet,  with  increased  turbulent  fluctuations  due  to  a  larger 
density  gradient,  ultimately  serving  to  prematurely  collapse  the  hot  jet.  The  energy  cascade  from  the  large 
bulk  mean  flow  to  the  smaller  scales  is  now  enhanced,  as  the  smallest  scales  are  now  relatively  larger  (as 
seen  in  the  structures  of  the  shear  layer  in  Figure  27). 

It  has  been  previously  shown  in  the  coupling  of  the  near  and  far-ficld  pressure,  that  the  column  mode 
exhibits  a  low  frequency  dependency,  which  is  prevalent  in  the  near-field  signature  of  the  hot  jet,  that  drives 
the  correlation  with  the  far-ficld  acoustics.  This  is  now  shown  to  be  the  case  in  the  source  velocity  field, 
as  well.  Results  presented  here  suggests  that  it  may  inherently  be  a  Reynolds  number  effect,  inhibiting  the 
jets  natural  development,  that  contributes  to  the  increased  low  frequency  noise  exhibited  in  the  hot  jet,  as 
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supported  by  Viswanathan.41 


H.  Noise  source  statistics 

The  turbulent  motion  of  compressible  fluid  is  responsible  for  the  generation  of  acroacoustic  noise.  In  the 
axisymmctric  jet,  this  noise  generation  is  governed  by  the  interaction  of  fluctuating  turbulent  quantities 
convecting  with  the  bulk  mean  flow.  Lighthill’s  wave  equation,  describing  the  relationship  between  the 
source  of  noise  in  free  shear  flows  and  the  resulting  distribution  of  external  fluctuations  can  be  solved  for 
the  fluctuating  density,  in  similar  fashion  as  in  Eq.l  described  previously. 


,  x  1  f  d2Tij  ,  . 


dy 


|x-y| 


(39) 


As  the  source  strength  in  the  acoustic  analogy  scales  as  a  double  divergence,  ,  it  follows  that  the  fluid 

motion  behaves  as  a  volume  distribution  of  acoustic  quadrupolcs.  Following  a  formulation  by  Proudman68 
and  Ribner,69  the  acoustic  power  output,  or  intensity,  in  the  far-ficld  defined  by, 


7(x)  = 


_2  _3 


^0  2 


Po 


(40) 


can  then  be  expressed  as  a  function  of  a  series  of  fourth-order,  quadrupole  space- time  correlations. 


7(x)  = 


XjXjXkXi 


a 


'  d2Tu  d2T, 


167r2aQpo^5  Jv  Jv  dr2  9t2 

where  only  the  perpendicular  longitudinal  axes  in  the  differentiation  are  significant,  scaling  with 


H^Mdydy' 

aJr,, 


(41) 


STa 

dxidxi 


Xj  1  dTjj 
x  a?,  dr2 


(42) 


The  ratio  Xi/x  is  the  direction  cosine  between  source  and  observer  location.  After  assumptions  neglecting 
viscous  and  temperature  effects,  i.e.  7);  =  puxUj,  the  resulting  acoustic  power  becomes, 


/(x)=TSSpXwwdy 


(43) 


Of  the  possible  thirty-six  components  contributing  to  the  quadrupole  correlations,  Ribner69  proposes  only 
nine  contribute  significantly  to  the  broadband  noise  pattern  of  the  jet.  Terms  involving  the  mean  velocity 
component  of  the  dominant  flow  direction  in  the  jet  deemed  shear  noise,  and  those  only  of  fluctuating 
quantities,  self  noise.  This  is  more  clearly  articulated  as  noise  contribution  from  turbulence- mean  flow 
interaction  and  turbulence-turbulence  interaction,  respectively.  The  contribution  of  each  to  the  distribution 
seen  in  the  far-field  intensity,  demonstrates  an  omnidirectional  pattern  (independent  of  $),  from  the  self 
noise  and  dipole-like  (Bcos2  0  4-  cos4  <£),  from  shear  noise.  The  combined  contribution  of  both  exhibits  an 
elongated  elliptical  pattern.  When  effects  of  convection  and  refraction  are  taken  into  account,  the  more 
common  ’figure-eight’  shaped  pattern  is  generated,  exhibiting  a  minimum  at  jet  centerline,  then  peaking  at 
shallow  aft  angle  and  diminishing  with  increased  aft  angle  separation.  The  rearward  angles  demonstrating 
a  minimum  intensity.  The  broadband  distribution  can  be  decomposed  to  evaluate  frequency  content  of  the 
these  noise  terms.  The  modeled  spectrums  reveal  that  the  self  noise  contribution  peaks  at  substantially 
higher  frequency  than  that  of  the  shear  noise,  resulting  from  a  quadratic  dependence  with  the  turbulent 
velocity  field  ( linear  relationship  with  shear  noise). 

Recalling  the  FS/LS  spectrum  of  Tam  et  al.23  it  is  apparent  that  the  two  theories  have  similar  founda¬ 
tions.  Tam  et  al.  suggesting  that  the  shape  of  the  far-field  spectrum  is  governed  by  the  dominant  large  scale 
contribution  at  shallow  aft  angles,  and  Ribner  demonstrating  that  the  largest  contribution  at  these  angles 
stems  from  shear  noise,  emanating  primarily  from  the  center  of  the  mixing  layer  (which  is  governed  by  large 
scale  motion,  driving  it  spreading). 
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The  possible  self-noise  and  shear-noise  terms  are  as  follows,  These  are  now  evaluated  to  asses  their 
contribution  in  both  the  cold  and  hot  jet.  Only  the  first  three  shear  noise  terms  are  considered,  as  the 
mean  flow  direction  (U\)  of  the  axisymmetric  jet  greatly  out  weighs  the  other  two  directions  ( radial  and 
azimuthal ),  giving  a  largely  more  significant  contribution.  Also,  following  the  reported  experimental  work 
of  Ukeiley  et  a/.,66  and  Bridges,60  only  the  normal  velocity  correlations  of  the  self-noise  terms  are  examined 
here,  as  these  have  been  shown  to  exhibit  the  largest  magnitudes  (u2u'2).  The  formulation  above  assumes 
isotropic  turbulence  and  neglects  the  contribution  of  the  second  source  term,  ( p  —  a^p),  which  has  been 
shown  to  scale  with  increases  in  temperature  fluctuation.  In  line  with  this  reasoning,  it  is  the  aim  of  this 
analysis  to  qualitatively  compare  the  contribution  of  the  Reynolds  stress  term(s)  in  both  jets  and  infer  as 
to  the  contribution  of  the  second  source  term  from  results.  The  overbar  in  the  above  terms  denotes  both 
a  spatial  and  temporal  average  of  two-point  space-time  correlations.  In  the  present  investigation,  the  time 
separation  is  fixed  at  r=0,  as  measurements  are  acquired  using  a  single  point  in  time,  multi-point  in  space 
PIV  system.  This  limits  the  ability  to  asses  the  frequency  relationship  of  individual  terms.  However,  it 
has  been  shown  in  the  numerical  analysis  of  Freund,71  that  maximum  correlation  are  exhibited  at  this  time 
separation. 

The  shear- noise  terms  representing  axial,  radial,  and  tangential  source  components  ( top  to  bottom),  at 
fixed  location  of  r/D= 0,  0.5,  &:  1  ( left  to  right),  at  a  streamwise  location  of  x/D= 4.5  are  seen  in  Figure  28 
for  both  jets.  As  seen  previously,  the  hot  jet  exhibits  a  larger  spatial  extent  in  nearly  all  contours.  Similar 
magnitudes  are  exhibited  in  both  jets,  with  a  maximum  along  the  centerline  and  decrease  seen  with  increasing 
outward  radial  distance.  This  is  consistent  with  results  reported  by  Jordan  &  Gervais,70  who  modeled  both 
isotropic  and  anisotropic  turbulence,  showing  that  the  shearing  of  turbulence  with  the  mean  flow  resulted 
in  longitudcnally  aligned  quadrupoles,  with  maxima  at  the  centerline  in  both  cases.  The  efficiency  of  these 
source  terms  depends  greatly  on  the  spatial  extent  over  which  the  turbulence  is  correlated.  It  has  been 
demonstrated  that  the  radial  distance  proves  to  be  the  most  significant  length  scale,  as  this  is  the  direction 
over  which  the  velocity  gradient  acts.  The  width  of  the  peaks  thus  signify  the  regions  of  most  efficient  sound 
generation.  As  such,  it  follows  that  the  hot  jet  would  exhibit  a  more  significant  contribution  to  the  far  field 
spectrum.  In  line  with  this  result,  the  hot  jet  contributes  more  low  frequency  acoustic  energy  in  the  direction 
of  mean  flow  (shallow  aft  angles). 

The  axial  source  components  are  shown  by  Jordan  &  Gervais70  and  Freund71  to  give  the  most  significant 
contribution,  however  Ukeiley  et  al.66  demonstrates  that  it  is  the  sources  associated  with  the  tangential 
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Figure  28.  Shear-noise  source  correlations  at  xjD— 4.5  ( normalized  by  f/^)(a)Cold  jet  (b)Hot  jet 
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Figure  29.  Self-noise  source  correlations  at  x/D= 4.5  ( normalized  by  U (a)Cold  jet  (b)Hot  jet 
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Figure  30.  Shear  noise  source  correlations  at  x/D= 8  ( normalized  by  U*t)  (a)Cold  jet  (b)Hot  jet 
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velocity  component,  143,  which  exhibits  a  maxima.  Here,  no  clear  trend  is  disccrnablc,  but  this  will  be 
revisited  following  the  presentation  of  the  remaining  sources.  The  self-noise  terms  offer  a  different  trend 
in  Figure  29,  as  the  cold  jet  exhibits  a  noticeably  larger  contours.  Here,  the  scaling  of  this  region  is  again 
significant  as  it  implies  an  acoustic  source  region,  however,  not  as  dependent  on  radial  distance  as  the  shear- 
noise  previously  exhibited.  These  sources  have  a  greater  dependency  on  turbulence  intensity.  Although  the 
magnitudes  arc  consistently  larger  in  the  hot  jet,  owing  primarily  to  the  increased  turbulent  fluctuations, 
the  region  of  effective  noise  generation  is  larger  in  the  cold  jet.  An  investigation  of  the  total  squared  energy 
under  the  curves,  as  intensity  scales  with  the  integrated  correlation,  reveal  a  higher  self-noise  contribution 
from  the  cold  jet.  This  suggests  that  the  cold  jet  is  a  more  efficient  high  frequency  noise  generator. 

Downstream  of  x/D=8,  the  contours  of  the  shear-noise  exhibit  similar  regions  over  which  the  turbulence 
is  correlated,  resulting  from  the  less  pronounced  differences  between  the  two  jets  here,  as  shown  previously. 
The  self-noise  sources  also  depict  a  similar  behavior,  with  a  more  pronounced  evidence  of  the  potential  core 
region  collapse  seen  in  the  hot  jet  (Figures  30  &:  31). 

With  downstream  separation,  at  r/D= 1,  the  self-noise  sources  suggest  a  breakdown  in  organization  of 
the  hot  jet’s  structure  with  downstream  evolution,  consistent  with  a  merging  of  the  mixing  layer.  A  similar 
trend  is  seen  in  the  cold  jet,  though  not  as  pronounced  ( structures  not  as  organized  overall  as  shown  by  POD 
analysis).  There  is  also  an  increase  in  the  magnitude  of  noise  contribution  and  spatial  extent,  in  both  jets, 
with  the  exception  of  the  centerline  where  ’peaks’  are  reduced  (overall  correlated  region  still  increases).  A 
seemingly  more  complicated  trend  is  seen  in  the  shear  noise  sources.  The  magnitudes  here  decrease  with 
downstream  distance  at  r/D= 0,  exhibit  a  slight  decrease  at  r/D= 0.5  ( more  consistent  in  the  hot  jet),  and 
increase  at  r/D=  1.  This  is  not  surprising  as,  these  terms  scale  with  the  local  mean  velocity.  A  review  of  the 
streamwise  profiles  in  Figure  16  highlights  the  reasoning  behind  the  observed  trend. 

In  each  jet,  at  both  streamwise  locations  observed,  the  self-noise  sources  contributions  outweigh  those 
of  the  shear- noise,  as  reported  by  others  71,70\66  There  is  a  less  pronounced  difference  exhibited  further 
downstream  of  the  potential  core.  Trends  regarding  the  dominance  of  individual  components  seen  elsewhere, 
is  not  clearly  discernablc  here.  Only  along  the  centerline  do  the  axial  components  (1,1)  consistently  outweigh 
the  radial  (2,2)  and  tangential  (3,3).  As  these  are  fourth  order  statistics  of  experimental  measurements,  the 
uncertainty  in  these  correlations  is  greatly  increased,  and  may  prove  to  give  vastly  different  results.  This  is 
especially  true  of  the  in-plane  ( radial  and  tangential)  components,  where  fluctuations  are  of  similar  order  as 
the  mean  value.  Acquiring  accurate  turbulence  statistics,  is  extremely  challenging.  However,  the  consistent 
trends  exhibited  in  the  values  reported  here,  incur  confidence  in  the  relative  behavior  of  the  source  terms. 
Thus  this  analysis,  in  the  least,  serves  as  a  good  qualitative  comparison  of  the  relative  contribution  of  the 
noise  sources  in  each  jet. 

As  the  contribution  of  the  second  source  term,  ( p  -  alp),  is  neglected  here  an  incomplete  picture  of  the 
source  strength  in  the  hot  jet  is  left.  However,  as  suggested  by  the  modeling  of  Morfey  et  a/.,40  and  shown 
numerically  by  Freund,71  it’s  contribution  serves  to  increase  the  radiation  of  noise  in  the  direction  of  shallow 
aft  angles.  ( This  is  even  apparent  at  the  low  Tr =0.86  in  Freund's  numerical  simulation).  In  this  regard,  the 
inclusion  of  this  term  would  serve  to  increase  the  already  dominant  low  frequency  shear-noise  (peak  radiation 
at  shallow  aft  angles)  of  the  hot  jet. 

These  results  are,  again,  consistent  with  the  findings  of  the  acoustic  far-field  analysis.  The  hot  jet  is  shown 
to  demonstrate  an  increased  contribution  of  low  frequency,  shear-noise,  which  is  predominantly  radiated 
in  the  direction  of  mean  flow.  The  far-field  spectrum  correspondingly  show  an  increased  contribution  of 
low  frequency  noise  at  shallow  aft  angles.  The  cold  jet,  conversely,  demonstrates  a  greater  contribution 
of  high  frequency,  self-noise,  at  a  similar  streamwise  location.  It  may  be  premature  here  to  suggest  that 
the  diminished  high  frequency  energy  in  the  acoustic  spectrum  of  the  hot  jet  is  tied  to  the  lesser  self-noise 
exhibited  without  further  analysis  over  a  larger  number  of  axial  locations,  however  the  correlation  is  striking. 
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Figure  31.  Self  noise  source  correlations  at  x/D=8  ( normalized  by  U^)  (a)Cold  jet  (b)Hot  jet. 
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I.  Conclusions 


The  acroacoustic  source  features  of  an  ’acoustically  matched’  Mach  0.6  jet,  with  temperature  ratio  Tr= 0.93 
and  1.7,  have  been  investigated.  With  exit  velocity  Uj  fixed,  and  core  flow  temperature  increased,  the  effect 
of  elevated  temperature  on  the  source  acoustics  was  isolated.  As  the  jet  Mach  number  ( Mj )  falls  to  0.44 
and  Reynolds  number  almost  half  of  the  cold  study,  2.8xl05  from  6.9xl05,  it  may  be  expected  that  less 
fluctuations  would  exist  as  the  jet  moves  towards  more  laminar  behavior.  However,  an  analysis  of  the  near- 
field  pressure  measured  by  an  azimuthal  array  positioned  near  the  jet  exit  (x/D= 2),  demonstrates  a  marked 
increase  in  low  frequency  energy,  seen  in  the  pressure  power  spectral  density  of  the  hot  jet,  accompanied 
by  a  slight  decrease  at  high  frequency.  A  modal  decomposition  reveals  the  presence  of  differing  dominant 
modal  contribution  in  each  jet,  ( column  mode  in  the  hot  jet  and  helical  in  the  cold).  The  more  distinct  peaks 
exhibited  in  the  coherency  curves  of  the  hot  jet,  suggests  a  more  consistent  scale.  Conversely,  broader  peaks 
in  the  cold  jet  imply  a  greater  distribution  of  scales  over  which  energy  is  correlated. 

The  mean  statistic  of  the  source  velocity  field  show  that  both  jets  have  similar  exit  velocity  and  a  similar 
spreading  rate  in  the  outer  shear  layer,  exhibiting  similar  half-widths  with  the  centerline.  However,  upon 
examination  of  the  shortened  potential  core  region  in  the  hot  jet  ( x/D=5  versus  6  in  the  cold ),  it  is  realized 
that  the  shear  layer  thickness  is  in  fact  larger,  owing  to  the  larger  collapsing  angle  (5.7°  in  the  hot  jet  and 
4.75°  in  the  cold).  The  turbulent  fluctuations  arc  seen  to  increase  in  all  velocity  components  in  the  hot  jet, 
with  a  predominance  in  the  axial  velocity  component,  at  nearly  3%.  The  data  from  both  jets  collapses  well 
using  r](x),  half  velocity  scaling,  as  well  as  the  Witze  correlation  parameter  for  jets  of  differing  temperature 
ratios,  which  gives  confidence  in  these  measurement. 

Reynolds  stresses,  arc  examined  in  both  jets,  again  revealing  trends  consistent  with  a  more  rapid  decay 
of  potential  core  region  of  the  hot  jet,  as  peaks  move  closer  toward  the  centerline  at  similar  downstream 
separation.  Examined  along  the  centerline  at  r/D— 0,  the  hot  jet  depicts  the  presence  of  larger  scales  closer 
to  the  jet  exit.  With  increased  downstream  distance,  the  differences  in  scaling  between  the  two  jets  exhibits 
only  subtle  differences,  owing  to  the  stunted  development  of  the  hot  jet  as  the  mixing  layers  have  merged,  no 
longer  drawing  energy  from  the  potential  core,  ( which  has  since  collapsed  further  upstream  in  the  hot  jet). 
The  length  scales  are  shown  to  be  consistently  larger  in  the  hot  jet,  as  compared  to  it’s  cold  counterpart,  with 
the  radial  extent  of  the  normalized  stresses  exhibiting  the  greatest  differences  in  all  terms,  an  effect  which  is 
shown  to  be  less  pronounced  with  increasing  distance  from  centerline.  The  tangential  normal  stresses,  serve 
as  an  exception  to  this  trend,  still  exhibiting  larger  scales  in  the  hot  jet,  but  in  azimuthal  variation. 

A  low-dimensional  analysis  of  the  two  source  velocity  fields,  using  the  ’snapshot’  POD  method,  reveals 
the  dominance  of  a  Fourier  column  mode-like  structure  in  the  hot  jet.  The  cold  jet  is  shown  to  exhibit  a 
dominant  Fourier  helical  mode-like  structure.  This  explains  the  observed  trend  seen  in  the  shear  stresses, 
whereby  the  column  mode  structure  of  the  core  exhibits  only  a  positive  lobe  with  the  helical  structure  of 
the  radial  component  ( helical  due  to  the  opposing  direction  of  motion  across  azimuthal  separation  due  to 
spreading  of  the  shear  layer).  Correspondingly,  the  turbulent  kinetic  energy  in  the  hot  jet  is  shown  to 
converge  more  rapidly,  depicting  a  higher  amount  of  the  system’s  total  energy  over  the  lower  ordered  modes, 
(more  pronounced  near  the  jet  exit).  The  largest  disparity  is  seen  at  x/D= 4.5,  where  the  first  POD  mode  of 
the  hot  jet  more  than  doubles  the  contribution  of  the  cold  (6%),  exhibiting  more  than  13%  of  the  flow  field’s 
total  energy.  A  rather  striking  result  is  observed  as  the,  POD  modal  structure  of  the  hot  jet  depicts  only 
clearly  well-organized  structures,  similar  to  those  of  Fourier  modes.  The  cold  jet  exhibits  some  deviation 
from  these  kinds  of  structures,  an  effect  which  is  show  to  exist  in  both  jets  at  downstream  location.  As  no 
periodicity  of  the  flow  field  is  presumed  in  the  decomposition,  this  trends  suggests  that  with  the  addition  of 
heat  the  flow  field  tends  to  a  more  radially  organized  structure,  resulting  in  a  more  homogeneous  and  periodic 
jet,  as  the  axisymmetric  column  mode  structure  is  shown  to  persist  longer  in  the  source  velocity  field  of  the 
hot  jet.  The  breakdown  of  this  organization  leads  to  the  trend  of  ’splitting’  Fourier-like  modes  observed, 
where  the  POD  modes  now  sum  to  form  Fourier  modes,  or  simply  exhibit  non-rcprcsentable  eigenfunctions. 
This  only  occurs  downstream  of  the  collapse  of  the  potential  core  in  the  hot  jet,  where  the  axisymmetric 
shear  layers  have  merged.  However,  at  fixed  acoustic  Mach  number,  this  may  well  be  a  lowered  Reynolds 
number  effect,  as  this  behavior  is  typical  of  a  more  ’laminar’  flow  structure. 

The  results  of  the  source  velocity  field  analysis,  are  in  line  with  those  reported  in  separate  analysis  of  the 
fluctuating  pressure  near-field.  The  increase  in  the  turbulence  levels  of  the  hot  jet  arc  sense  by  the  near-field 
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array  suggesting  that  the  region  surrounding  the  flow  has  a  strong  dependency  34 1  53 \35  The  dominant 
column  mode  of  the  hot  jet  is  sensed  as  an  increased  low  frequency  contribution  in  the  pressure  near-field. 
The  distinct  peaks  of  the  pressure  coherency  spectrum,  suggesting  the  presence  of  a  more  consist  scale  in  the 
hot  jet  and  larger  distribution  of  scales  in  the  cold,  is  verified  through  the  low-dimensional  reconstruction  of 
the  dominant  structures  of  the  source  velocity  field.  The  more  gradual  change  in  the  ’richness’  (shapes  and 
magnitudes  of  the  peaks),  of  the  flow  features  seen  in  the  cold  jet,  with  the  addition  of  a  similar  number  of 
POD  modes,  implies  a  similar  weighting  of  modes.  In  the  hot  jet  the  energy  distribution  is  such  that  once 
the  high  energy,  bulk  flow  features  arc  realized,  the  higher  ordered  modes  contribute  little  to  the  shape  or 
magnitude  of  the  curves,  until  at  a  sufficiently  large  mode  count. 

The  acoustic  features  of  the  each  jet  as  sensed  by  the  far-field  array  of  microphones  exhibits  some 
distinction,  in  line  with  findings  of  others38,60,27’.56  The  directional  dependence  of  the  far-field  spectrum 
with  frequency  is  highlighted,  as  the  spectral  densities  reveal  an  increase  in  the  amplitude  of  the  hot  jet 
at  nearly  all  observer  angles,  with  a  marked  increase  in  the  low  frequency  energy  at  shallow  aft  angles. 
The  overall  sound  pressure  levels  demonstrate  a  2 dB  increase  seen  at  the  shallowest  aft  angle  of  0=15° 
and  a  slight  decrease  of  nearly  a  half  dB  seen  at  the  <p=90°  location.  This  increase  is  a  result  of  increased 
low  frequency  noise  contribution  in  the  hot  jet,  coupled  with  a  decrease  at  high  frequency,  as  compared 
to  the  cold  jet,  seen  in  the  acoustic  spectral  densities.  An  analysis  of  the  frequency  contribution  of  each 
microphone  location  reveals  more  pronounced  changes  in  amplitude  with  increasing  frequency  at  shallow 
aft  angles,  substantiating  the  trend  of  low  frequency  noise  dominance  at  shallow  aft  angles,  reported  in 
the  OASPL  increase.  Correspondingly  the  reduction  in  high  frequency  noise,  is  reflected  in  the  diminished 
contribution  in  OASPL  of  the  microphone  located  at  the  largest  observer  angle  of  0=90°. 

The  far-field  acoustic  features  arc  shown  here  to  be  reflective  of  the  changes  in  the  source  velocity  field, 
similarly  to  the  near-field  pressure.  The  marked  increase  in  low  frequency  noise  at  shallow  aft  angles,  can 
be  linked  to  the  dominant  low  frequency  column  mode  present  in  the  hot  jet.  In  this  similar  regard,  the 
reduction  of  the  higher  frequency  noise  can  be  linked  to  the  increase  in  the  scales  of  the  structure  ( small  scales 
producing  high  frequency  noise).  Thus,  the  shift  to  a  dominant  low  frequency  column  mode  observed  in  both 
the  near-field  pressure  signature,  along  with  the  increase  in  low  frequency  acoustic  energy  seen  predominantly 
at  shallow  aft  angles,  again,  serves  as  evidence  of  the  the  large  scale’s  contribution  predominantly  to  the 
shallow  aft  angles. 

The  increase  in  OASPL  demonstrates  a  greater  influence  on  low  frequency  energy  increase,  than  high 
frequency  reduction.  Where  noise  production  is  concerned,  the  Reynolds  stresses,  viewed  as  a  simplified 
representation  of  the  Lighthill’s  stress  tensor,  Tij  =  pUiUj ,  exhibit  increased  magnitudes  in  the  hot  jet.  As 
the  turbulence  levels  do  not  show  significant  differences,  the  link  between  these  stresses  and  the  increase 
in  OASPL  is  weakened  by  the  vast  difference  in  mean  density  of  each  jet  (pCoid= 1-428  and  phot=0.759). 
Therefore,  the  Reynolds  stress  source  term’s  contribution  to  noise  generation  in  the  hot  jet  would  no  doubt 
be  diminished. 

The  mechanism  responsible  for  the  additional  low  frequency  noise  is  sought  in  an  investigation  of  the 
acoustic  sources,  using  the  quadrupolc  correlations.  These  fourth  order  space-time  statistics  are  calculated 
directly  from  the  fluctuating  source  quantities.  The  temporal  domain  is  suppressed,  at  t— 0,  due  to  the  use 
of  PIV  to  acquire  the  velocity  field.  The  shear  noise  terms,  contributing  to  low  frequency  fluctuations,  are 
shown  to  depict  similar  magnitude  in  each  jet,  with  a  maximum  along  the  centerline.  As  it  is  the  radial 
extent  of  the  correlation  which  greatly  influences  the  efficiency  of  these  source  terms  ( velocity  gradient  acts 
over  spatial  extent),  the  width  of  the  peaks  thus  signify  the  regions  of  most  efficient  sound  generation.  The 
hot  jet,  with  larger  sealed  structures  exhibiting  a  greater  correlation  length  scale,  in-turn  demonstrates  a 
more  significant  contribution  to  the  far  field  spectrum.  The  peak  angle  of  radiation  as  modeled  by  Ribner  69 
for  these  terms  is  in  the  direction  of  mean  mean  flow,  or  at  shallow  aft  angles. 

The  self  noise  terms  show  a  marked  contribution  in  peak  intensity  in  the  cold  jet,  exhibiting  noticeably 
larger  spatial  extent  of  correlation.  As  these  sources  have  a  greater  dependency  on  turbulence  intensity,  the 
magnitudes  arc  consistently  larger  in  the  hot  jet  ( owing  primarily  to  the  increased  turbulent  fluctuations). 
However,  intensity  scales  with  the  integrated  correlation  over  the  region  of  effective  noise  generation,  which  is 
larger  in  the  cold  jet,  revealing  a  higher  self-noise  contribution  from  the  cold  jet.  The  cold  jet  is  thus  a  more 
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efficient  high  frequency  noise  generator.  These  terms  arc  know  to  exhibit  an  omnidirectional  dependence.  At 
each  measured  location,  the  self-noise  sources  contribution  consistently  outweighs  those  of  the  shear-noise. 
There  is  a  less  pronounced  difference  exhibited  further  downstream  of  the  potential  core.  However,  as  the 
the  self  noise  terms  illustrate  the  contribution  of  high  frequency  energy  to  the  acoustic  spectrum,  then  the 
contribution  of  theses  terms  must  outweigh  the  low  frequency  noise  of  the  shear  noise  source  terms,  in  both 
jets. 

In  a  similar  line  of  reasoning,  the  diminished  contribution  of  high  frequency  noise  in  the  hot  jet  should  be 
less  pronounced  than  the  increase  in  low  frequency  noise,  as  the  acoustic  spectra  demonstrates  a  larger  low 
frequency  change.  As  no  clear  distinction  supporting  this  can  drawn  from  the  quadrupole  correlations  (shear 
noise  sources  in  both  jets  depicting  similar  magnitudes  with  only  a  slightly  larger  spatial  extent  in  the  hot),  the 
contribution  of  the  entropy  source  term  is  deemed  a  significant  contributor.  The  formulation  used  to  derive 
these  terms  does  not  include  the  contribution  of  the  second  source  term  (p  —  a2p)6ij,  however  as  evaluated 
numerically  by  Freund,71  it’s  contribution  is  shown  to  exhibit  a  marked  increase  in  noise,  predominantly  at 
shallow  aft  angles.  This  result  is  found  in  the  cold  jet  at  Tr=0.86.  As  such,  it  may  be  inferred  that  the 
contribution  of  this  term  will  be  more  pronounced  in  the  case  where  a  temperature  gradient  exists,  (higher 
temperature  ratio).  This  contribution  would  serve  to  increase  the  already  dominant  low  frequency  shear-noise 
contribution  (in  the  direction  of  peak  radiation  at  shallow  aft  angles)  of  the  hot  jet,  as  compared  to  the  cold. 


The  presence  of  the  ’additional  source’  of  acoustic  energy  in  the  hot  jet  can  now  be  described  as  resulting 
from  the  increased  organization  of  the  velocity  field,  causing  the  column  mode  structure  to  persist  longer  than 
its  cold  counterpart.  Larger  scales  are  prevalent,  owing  to  the  stagnation  of  the  natural  development  of  the 
jet,  and  the  increased  turbulent  fluctuations  (driven  by  a  larger  density  gradient  between  core  and  entrained 
fluid).  This  serves  to  prematurely  collapse  the  hot  jet.  The  scaling  of  these  structures  contributes  greatly 
to  the  far-ficld  spectrum  in  the  form  of  low  frequency  shear  noise,  radiating  primarily  in  the  direction  of 
propagation,  at  shallow  aft  angles.  The  resultant  field  is  loosely  described  as  a  more  laminar  jet  with  increased 
turbulent  fluctuations.  This  suggests  that  it  may  inherently  be  a  Reynolds  number  effect ,  as  this  trend  of 
’laminarization’  is  associated  with  a  decrease  in  the  effective  Reynolds  number  in  the  fixed  Mach  number 
experiment,  as  supported  by  Viswanathan.41  (Jet  Mach  number  (Mj)  falls  to  0.44  form  0.6  and  Reynolds 
number  almost  half  of  the  cold  study,  2.8xl(P  from  6.9xl(P).  However,  the  less  pronounced  difference  in  the 
source  contribution  of  each,  docs  speak  to  the  increased  contribution  of  the  second,  entropy  source  term. 

1 .  Discussion 

This  study  comprises  several  independent  analysis,  i.e.  the  source  velocity  field,  fluctuating  pressure  sur¬ 
rounding  it,  and  the  acoustic  far-field  response.  It  has  been  demonstrated  that  the  fluctuating  pressure, 
sampled  here  at  x/D= 2,  r/D  ~1,  accurately  captures  the  characteristic  dynamics  of  the  source  velocity 
field34,53’.35  In  turn,  it  is  also  shown  that  the  acoustic  far-field  exhibits  a  strong  relationship  with  the 
fluctuating  near-field  signature.52 

Next  steps  would  include  coupling  the  source  velocity  field  with  the  surrounding  near-field  pressure,  via 
the  application  of  the  modified  complementary  technique,  as  detailed  by  Bonnet  et  al.42  As  the  technique 
relies  greatly  on  the  measure  of  correlation  between  the  two  fields,  the  current  results,  establishing  a  re¬ 
lationship  between  the  two,  are  promising.  Once  a  reliable  set  of  estimation  coefficients  are  obtained,  the 
undcrsamplcd  PIV  velocity  field  can  be  estimated  to  realize  a  time- resolved  field.  This  in  turn  provides  a 
means  by  which  the  simplified  stress  term  in  Lighthill’s  acoustic  analogy,  Tij  =  puxUj ,  cam  be  calculated 
and  ultimately  used  to  estimate  the  far-field  acoustics.  With  the  sampled  acoustic  field  already  available, 
comparisons  can  be  made  between  the  two  to  determine,  firstly,  whether  the  assumption  of  neglecting  the 
entropy  term  (p  -  aQp)Sij  is  justifiable  in  the  cold  jet.  Secondly,  this  term’s  contribution  in  the  hot  jet  can 
be  evaluated,  as  the  sampled  acoustic  far-ficld  of  the  cold  jet  could  serve  as  a  benchmark,  whereby  levels 
above  this  would  be  deemed  additional  noise. 

This  analysis  highlights  the  amount  of  information  that  is  lacking  as  a  result  of  the  inability  to  sample 
inflow  pressure  and  density  directly.  Panda  et  al.72  has  recently  applied  a  Rayleigh  scattering  technique  to 
get  a  measure  of  the  density  fluctuations,  however  with  current  technology,  measurements  of  inflow  pressure 
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would  require  bulky  probes  resulting  in  the  contamination  of  the  ’true’  acoustic  far-ficld  signature.  The  use 
of  numerical  simulations  such  as  Large  Eddy  Simulation  and  Direct  Numerical  Simulation  provide  solutions 
to  evaluating  the  contribution  of  the  entropy  term.  These,  however  rely  greatly  on  the  initial  and  boundary 
conditions  employed  to  simulate  the  flow  field.  As  well,  extending  the  computational  domain  to  the  far-ficld 
requires  further  assumptions  to  be  made. 

Where  control  of  noise  is  concerned,  these  results  indicate  that  it  is  the  intrinsic  structure  of  the  jet  that 
must  be  manipulated.  The  increased  low  frequency  noise  is  seen  to  result  from  differences  associated  with 
the  scales  and  modal  ordering  of  the  dominant  structures  in  each  jet.  A  delicate  balance  must  be  reached 
in  manipulating  these  through  the  use  of  upstream  perturbation  or  fluid  injection.  Column  mode  forcing  at 
the  jet’s  preferred  instability  serves  to  increase  the  low  frequency  noise,  whereas  at  higher  modes,  a  small 
reduction  can  be  achieved  as  a  result  of  the  combined  effect  of  low  frequency  reduction  and  high  frequency 
increase,  Samimy.75  It  is  in  this  regard  that  active  flow  control  schemes  must  be  adapted  to  include  more 
complex  forcing  signals  or  strategics  which  couple  differing  modes,  in  hopes  of  achieving  a  more  broadband 
reduction  in  noise. 


V.  Dynamical  Systems  Modeling 

Several  attempts  have  been  made  at  investigating  the  dynamics  of  the  high-speed  jet,  both  quantitatively 
and  in  a  time-resolved  manner.  This  involves  measuring  flow  quantities  so  as  to  capture  from  the  smallest  to 
the  largest  time-scale  events  in  the  flow  field.  Prom  these  type  of  experimental  or  computational  measure¬ 
ments,  dynamical  information  can  be  extracted  from  which  much  can  be  learned  as  to  the  types  of  coherent 
processes  that  occur  and  how  they  relate  to  the  far-field  sound. 

Since  the  1970’s  it  is  widely  recognized  that  even  in  very  high  Reynolds  number  flows,  there  exists 
coherent  structures  that  demonstrate  distinctive  features  such  as  organized  vorticity  and  exhibit  repeating 
patterns  of  turbulent  motion4,16’.15  Unlike  previous  beliefs  that  turbulent  flows  were  dominated  by  random 
motion  and  chaos,  through  advanced  measurements,  simulations,  and  analytical  data  processing  tools,  such 
features  arc  now  recognized  in  highly  turbulent  flows  as  well.  The  problem  lies  in  discerning  a  physical 
interpretation  of  these  results  to  explain  the  mechanisms  at  work  and  build  models  to  represent  them.  Such 
models  arc  important  in  furthering  our  understanding  of  turbulent  motion,  as  numerically  solving  the  full 
Navier-Stokes  equations  for  high  Reynolds  number  flows  is  not  a  realizable  goal  in  the  foreseeable  future 
of  computing  evolution.  The  development  of  dynamical  systems  models  can  therefore  improve  on  Reynolds 
Averaged  Navier-Stokes  (RANS)  techniques  or  validate  Large  Eddy  Simulations  (LES)  of  such  flows. 

Direct  Numerical  Simulations  (DNS)  of  high-speed  jets  performed  at  very  low  Reynolds  number,  along 
with  computations  of  the  instantaneous  aeroacoustic  field,  have  provided  great  insight  into  the  dynamics 
of  these  flows  and  their  relation  to  the  far-ficld  sound.  Freund  et  a/.11  have  applied  flow  control  using 
adjoint  methods  in  a  DNS  simulation  of  a  two-dimensional  high-speed  jet  with  the  control  objective  of 
minimizing  the  far-field  sound.  The  control  is  performed  in  a  ’black-box’  near  the  jet  exit  and  modifications 
to  the  control  input  arc  made  iteratively  at  retarded  times  to  slightly  modify  the  dynamics  of  the  flow, 
minimizing  the  sound  propagated  to  the  far-field.  A  quiet  jet  was  in  this  manner  simulated  and  compared 
to  the  uncontrolled  loud  jet.  Quite  surprisingly,  a  simple  comparison  of  the  instantaneous  dynamics  of  the 
controlled  to  the  uncontrolled  flow  does  not  show  a  significant  difference  in  the  organization  of  the  large-scale 
structures  of  the  flow,  which  indicates  that  the  noise  production  mechanism  is  quite  subtle.  Only  once  a 
processing  of  the  data  using  proper  orthogonal  decomposition  (POD)  is  performed  does  a  difference  in  spatial 
phase  appear  between  the  first  modes.  This  type  of  study  of  the  dynamics  of  jets  shows  that  a  similarly 
subtle  control  strategy  needs  to  found,  when  such  a  slight  change  in  the  flow  has  such  a  great  impact  on  the 
far-field  sound. 

Early  experimental  measurements  of  the  dynamics  of  high-speed  jets  were  performed  using  rakes  of 
hot-wires.  Multi-point  correlations  of  velocity  components  coupled  with  assumptions  of  homogeneity  and 
periodicity  in  the  jet  flow  flied,  allowed  the  application  of  low-dimensional  techniques  such  as  Stochastic 
Estimation  and  POD.1  Citriniti7  took  measurements  of  the  dynamics  of  the  axisymmetric  incompressible 
jet  at  one  downstream  position  using  an  in-house  designed  rake  of  138  hot-wires.  The  experiment  was 
then  carried  out  at  multiple  near-field  positions  by  Gamard26  and  multiple  far-ficld  positions  by  Jung,7  to 
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study  the  evolution  of  the  azimuthal  modes  as  the  jet  develops  downstream.  Computational  efforts  towards 
measuring  the  dynamics  of  axisymmetric  jet  flows  arc  today  feasible  at  rather  low  Reynolds  numbers  (5xl03) 
in  direct  simulations71  and  at  intermediate  Reynolds  numbers  (5xl04)  using  LES.56  These  investigations  have 
the  advantage  of  containing  time-resolved  information  of  multiple  flow  quantities  simultaneously  (velocity, 
density,  acoustics,  hydrodynamics,  vorticity,  etc)  from  which  much  valuable  information  can  be  gained. 
However  these  methods  arc  limited  by  the  rather  low  Reynolds  numbers  attainable  and  the  small  amount  of 
time  scries  computed. 


A.  Measuring  flow  dynamics 

Experimental  time- resolved  measurements  of  turbulent  flow  fields  certainly  emerged  with  the  advent  of  hot¬ 
wire  anemometry,  principles  of  which  date  back  to  the  beginning  of  the  twentieth  century.  This  proved  to 
be  an  extremely  effective  and  accurate  tool,  but  limiting  in  terms  of  application.  At  higher  velocities,  the 
intrusive  nature  of  the  technique  would  cause  probes  to  break  under  the  high  dynamic  pressure  loading, 
and  generate  turbulent  vibrations.  With  the  introduction  of  lasers  came  laser  Doppler  anemometry  (LDA), 
which  although  technically  much  more  challenging  and  costly  for  the  experimenter,  presented  the  advantage 
of  being  a  non-intrusivc  technique.  It  was  now  possible  to  obtain  3-component  single-point  measurements 
of  velocity  in  high-speed  flows,  without  disturbing  the  flow.  As  this  optics  based  measurement  is  based  on 
the  Doppler  shift  provoked  by  a  particle  crossing  the  measurement  volume  scattering  light  with  a  frequency 
shift  proportional  to  its  velocity,  its  main  limitation  is  flow  seeding.  Non-intrusive  multi-point  velocity 
measurements  were  only  developed  toward  the  late  1980s  when  high  powered  lasers  and  fast  response  cameras 
were  combined  to  form  particle  image  velocimetry  (PIV).  This  technique  enabled  accurate,  spatially  well- 
resolved  velocity  measurements  simultaneously  in  an  entire  plane.  Again,  however,  this  technique  proved  to 
also  be  limited,  as  temporal  resolution  in  high-speed  flow  was  unattainable.  However,  much  can  be  gained 
from  such  measurements  if  one  has  the  desire  to  study  the  statistical  properties  of  a  flow  such  as  mean 
velocity  profiles,  phase- averaged  velocity  fields,  or  simply  collect  databases  of  instantaneous  snapshots  of 
velocity  fields  to  then  perform  multi-point  spatial  correlations  and  other  post-processing  techniques. 

Experimentally  measuring  the  dynamics  of  the  high-speed  jet  flow  can  prove  a  challenging  task.  We 
arc  today  experimentally  limited  to  two-dimensional  plane  measurements  (3  components).  However,  planar 
measurements  can  be  very  informative,  especially  when  all  three  components  of  velocity  can  be  measured 
simultaneously  on  a  relatively  fine  grid,  such  as  when  using  stereoscopic  particle  image  velocimetry  (PIV). 
This  optics-based  measurement  technique  is  attractive  not  only  because  of  its  non-intrusivcncss  (the  flow 
remaining  undisturbed  by  the  measurement)  but  also  due  to  its  range  of  application.  It  is  not  prudent  to  use 
rakes  of  hot-wires  in  high-speed  compressible  flows.  PIV  is  therefore  the  only  experimental  instrument  capa¬ 
ble  of  providing  planar  velocity  data  at  fine  spatial  resolution  in  high-speed  flows  such  as  the  axisymmetric 
jet.  The  next  step  in  dynamical  measurements  of  such  flows  was  the  recent  development  of  the  multi-plane 
PIV  measurement  technique7  which  presents  the  advantage  of  being  a  non-intrusivc  technique  applicable  to 
many  types  of  flows.  This  enables  a  direct  measurement  of  the  three-component  acceleration  field.  Depend¬ 
ing  on  application,  Eulcrian  or  Lagrangian  acceleration  can  be  realized  by  finite  difference  calculation  of  the 
consecutive  PIV  measurements;  the  latter  being  delayed  by  a  carefully  chosen  time  lapse  that  is  tied  to  the 
time  scales  of  the  flow  field  under  consideration. 


B.  Dual-time  PIV  experiment 

The  pioneering  use  of  multiple  PIV  systems  was  proposed  and  later  implemented  by  Kahlcr  &:  Kompcnhans' 
to  develop  a  tool  capable  of  measuring  time  derivatives  of  the  velocity  (z.e.  acceleration)  as  well  as  spatial 
derivatives  (z.e.  velocity  gradient  tensors,  strain  rate  tensors  and  vorticity  field)  in  a  whole  plane.  Several 
research  groups  have  since  performed  dual-plane  PIV  using  two  PIV  systems  for  measurements  in  two  parallel 
and  adjacent  planes;  as  well  as  dual-time  PIV  using  both  systems  in  the  same  measurement  plane  separated 
in  time  by  a  short  delay.  The  latter  technique  enables  measurement  of  acceleration  or  more  generally  of 
the  time  derivative  of  the  velocity  field  which,  associated  with  the  proper  orthogonal  decomposition  (POD), 
enables  dynamical  systems  development  with  a  potential  use  for  closed-loop  flow  control. 
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In  the  present  configuration  where  the  consecutive  PIV  measurements  arc  performed  at  the  same  location, 
the  Eulerian  acceleration  (or  local  acceleration)  is  computed  directly.  The  Lagrangian  acceleration  (or 
material  acceleration)  is  the  whole  left-hand  side  of  the  Navicr-Stokcs  equations  defined  as, 


Du 

~Dt 


du 

di 


+  (u  •  V)u 


(44) 


the  first  term  on  the  right-hand  side  being  the  Eulerian  acceleration  and  the  second  term  being  the  non-linear 
convective  acceleration.  To  measure  this  quantity  experimentally,  one  has  to  follow  each  particle  as  it  moves 
through  the  velocity  field  to  take  into  account  both  contributions  to  the  acceleration  of  the  particle,  the 
unsteadiness  of  the  velocity  field  and  the  change  in  velocity  as  the  particle  moves  through  this  field.  Two 
stereoscopic  PIV  systems  arc  used  in  the  present  experiment  in  an  identical  cross-flow  plane  perpendicular 
to  the  jet  axis,  capturing  all  three  components  of  acceleration.  Using  finite  differencing,  local  acceleration  is 
computed  for  each  pair  of  velocity  field  obtained.  The  process  is  reiterated  for  different  downstream  positions 
ranging  from  3  to  10  jet  diameters  so  as  to  obtain  a  3-dimensional  acceleration  field  of  the  high-speed  jet. 
A  picture  of  the  setup  is  shown  in  Figurcrcffigidual. 


Figure  32.  Dual-time  PIV  experiment  setup. 


The  camera  arrangement  purposely  makes  the  angle  between  the  cameras  in  each  pair  equal  to  make 
the  measurements  from  both  systems  as  identical  as  possible.  This  minimizes  the  differences  in  the  level 
of  error  which  would  be  amplified  when  computing  the  time  derivatives.  The  light  sheets  are  orthogonally 
polarized  and  linear  polarizing  lenses  arc  mounted  on  each  camera  to  ensure  separate  measurements.  Special 
care  and  time  was  taken  in  making  the  four  laser  sheets  identical  and  aligning  them  to  minimize  the  error 
between  the  measurements  of  both  systems.  The  sheet  thickness  is  3  mm  which  is  thick  enough  to  enable 
the  measurement  of  the  out-of-plane  component  of  velocity  (uz),  which  requires  that  a  particle  traveling 
through  the  light  sheet  is  illuminated  by  both  pulses  of  the  lasers.  In  the  present  ease,  the  delay  between 
pulses  is  set  to  5  /xs,  which  allows  a  particle  with  a  velocity  of  200  m/s  to  move  1  mm  within  the  3  mm 
thick  laser  plane. 

The  choice  of  A t  is  crucial  in  the  finite  difference  computation  of  the  time  derivative  of  the  velocity.  If 
A t  is  too  small,  An  =  ix(z,  t  +  At)  -  u(x,  t)  will  be  of  the  same  order  of  magnitude  as  the  error  associated 
with  the  PIV  measurements  themselves.  If  At  is  too  large,  then  the  dynamics  of  the  flow  at  the  time  scales 
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of  interest  can  not  be  resolved.  Based  on  an  extensive  study  on  this  subject  was  carried  out  a  priori,  a  At 
of  25  /is  was  chosen  for  the  downstream  positions  from  x/D= 3  to  x/D= 8  and  a  Au  of  35  /rs  for  x/D=S. 5 
to  x/D=  10,  as  the  time  scales  in  the  flow  arc  larger  further  downstream. 

C.  Low-order  dynamical  system  development 

Given  a  complex  system  fluid  flow,  the  goal  in  developing  reduced  order  models  (ROM)  is  to  solve  only  a  part 
of  the  infinite  dimensional  governing  evolution  equations  (partial  differential  equations)  by  projecting  these 
equations  onto  a  subspacc  of  finite  dimensions,  made  up  of  a  reduced  basis  of  empirical  eigenfunctions  (like 
those  provided  by  the  POD).  The  level  of  complexity  of  a  fluid  flow  increases  with  Reynolds  number,  but 
experiments  at  many  different  Reynolds  numbers  have  shown  that  low-order  analysis  is  effective  throughout 
the  whole  range  and  that  organized  structures  exist. 

A  Galcrkin  projection  is  the  projection  of  a  system  of  partial  differential  equations  onto  a  finite  basis  of 
linearly  independent  functions.  This  projection  can  then  be  truncated  to  model  only  a  part  of  the  system. 
The  proper  orthogonal  decomposition  is  a  method  that  conveniently  complements  the  Galcrkin  projection 
since  the  POD  basis  is  orthonormal  by  construction  and  the  functions  are  optimized  in  terms  of  the  square 
of  the  velocity  (i.c.  the  energy),  hence  when  truncating  the  system  at  a  low  mode  number,  only  the  most 
energetic  part  of  the  flow  is  modeled  in  the  system.  Since  an  empirical  method  is  used  in  the  present  work 
and  the  full  Galcrkin  projection  with  explicit  calculations  of  the  projection  coefficients  is  not  performed,  the 
analysis  tied  to  the  Galcrkin  projection  in  the  case  of  incompressible  flows  is  summarized  in  Bergmann.7  The 
equations  resulting  from  the  projection  are  quadratic  (as  are  the  Navier-Stokes  equations)  but  an  important 
point  lies  in  the  fact  that  the  Reynolds  decomposition  is  used  in  the  present  case.  Indeed,  when  separating 
the  mean  flow  from  the  fluctuations  in  the  projections,  extra  terms  describing  the  interaction  between  the 
mean  velocity  and  the  turbulent  fluctuations  arc  derived  by  Ukeiley  et  a/.14  Since  these  terms  arc  cubic 
in  nature,  an  additional  cubic  term  will  be  added  to  the  modeling  ODE  to  model  this  interaction.  The 
second  motivation  for  adding  a  cubic  term  in  our  case  is  that  the  flow  cannot  be  considered  incompressible, 
even  though  the  previous  developments  were  made  using  the  incompressible  version  of  the  Navier-Stokes 
equations.  Rowley  et  al.8  apply  a  Galcrkin  projection  to  the  full  compressible  Navier-Stokes  equations 
and  introduce  a  more  appropriate  inner  product  adapted  for  compressible  flows.  They  also  show  that  since 
the  derived  compressible  equations  are  this  time  cubic  in  nature,  a  cubic  model  needs  to  be  considered. 
Therefore,  the  following  general  ODE  will  be  used  in  the  development  of  an  empirical  low  order  dynamical 
system  of  the  jet,  as  was  similarly  used  by  Perret,21 

^  NloDS  NlODS  NlODS  NloDS  NLODS  NlODS 

—J—  =  Dn  -f  )>  ]  Lnjdj  +  }  '  ^  "  Qnjkajak  +  ^  ^  ^  ]  CnjklQ>jQ>kal  (45) 

at  j=  1  j=  1  k= 1  j  =  1  k=  1  /= 1 

where  Dn,  Lnj ,  Qnjk  and  Cnjki  are  the  constant,  linear,  quadratic  and  cubic  terms  of  the  dynamical  system 

that  need  to  be  identified.  As  discussed  above,  these  can  be  calculated  directly  using  the  analytical  expres¬ 
sions  derived  from  the  Galerkin  projection.  Very  finely  resolved  data  (both  in  time  and  space)  is  however 
needed  for  such  calculations  and  this  direct  technique  would  not  be  appropriate  in  the  case  of  experimental 
PIV  data  where  a  significant  amount  of  noise  would  be  introduced  by  the  double  spatial  derivatives  and 
the  computation  of  higher  order  terms.  These  coefficients  will  here  be  identified  using  an  empirical  method 
commonly  referred  to  as  the  ’moments’  method  developed  at  the  Laboratoirc  d’Etudes  Aerodynamiques  at 
the  Univerity  of  Poitiers,  France.  The  present  work  extends  this  technique  and  its  application  to  a  much 
higher  speed  flow  to  see  how  a  more  turbulent  and  less  ’coherent’  database  affects  the  prediction  capability 
of  this  method. 

The  moments  method  only  requires  an  ensemble  of  statistically  independent  samples  of  the  measurement 
and  not  a  time- resolved  signal.  The  method  in  fact  converges  faster  when  trained  with  independent  samples 
rather  than  with  slightly  correlated  measurements.  The  input  database  comprises  for  each  independent  time 
tx  a  pair  (an(ti)ydan/dt(ti))\  i  =  1,  2,  ...,  N.  A  linear  forward  differencing  algorithm  is  used  to  compute  the 
time  derivative  of  the  POD  coefficients, 

do>n  /,  \  &n(^0  At)  0,n(t o)  / 

W(to)  = - At -  (46) 
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The  advantage  of  using  DT-PIV  is  that  this  term  is  measured  directly.  To  determine  each  one  of  the  con¬ 
stant,  linear,  quadratic  and  cubic  coefficients  one  needs  to  derive  as  many  equations  as  there  arc  unknowns. 
In  this  aim,  several  moments  of  Eq.  45  are  taken  and  ensemble  averages  are  performed  using  the  experi¬ 
mental  database.  The  size  of  the  problem  can  become  very  quickly  unpractical  especially  when  including 
cubic  terms.  One  of  the  problems  inherent  to  low-order  dynamical  systems  in  fluid  dynamics  is  that  the 
truncation  performed,  greatly  reduces  the  dimensions  of  the  problem  but  also  neglects  most  of  the  smaller 
scales  responsible  for  dissipating  the  energy  through  viscous  forces.  A  common  resulting  behavior  is  for  the 
dynamical  system  to  rapidly  diverge.  To  model  the  truncated  scales  and  drain  the  energy  in  the  system,  an 
ad  hoc  artificial  viscosity  is  introduced  to  consistently  adjust  the  linear  coefficient  throughout  all  the  modes. 
The  number  of  modes  kept  in  the  low-order  system  is  therefore  always  chosen  to  ensure  the  problem  is  wrcll 
conditioned.  This  issue  restricted  the  number  of  modes  retained  in  the  modeling  to  Nlods  =  65  in  the 
quadratic  ease  and  Nlods  =  12  in  the  full  modeling  with  the  cubic  terms. 

A  total  of  6  low-order  dynamical  systems  arc  developed  at  the  downstream  position  of  x/D= 8  to  study 
the  influence  of  the  different  modeling  parameters  on  the  prediction  of  the  dynamics  and  hence  to  find  which 
are  critical  to  building  a  successful  dynamical  system  in  a  highly  turbulent  flow.  These  arc  detailed  in  Table6 
below. 


Table  6.  Parameters  used  in  the  6  different  low-order  dynamical  systems. 


System 

Order 

Nlods 

N. 

*0 

Total  Eqns. 

%  TKE 

Predicted  time 

LODS  1 

cubic 

12 

4000 

23 

455 

30 

10ms 

LODS  2 

cubic 

12 

4000 

35 

455 

30 

10ms 

LODS  3 

quadratic 

30 

749 

23 

496 

47 

25ms 

LODS  4 

quadratic 

48 

3000 

23 

1225 

57 

10ms 

LODS  5 

quadratic 

65 

4000 

23 

2211 

64 

10ms 

LODS  6 

quadratic 

65 

4000 

25 

2211 

64 

10ms 

The  first  two  LODS  were  developed  from  the  full  cubic  ODE  and  from  a  modeling  stand  point,  contains 
the  higher  order  terms  necessary  to  model  the  interactions  between  the  mean  flow  and  the  fluctuating 
velocities,  since  the  Reynolds  decomposition  is  used.  Both  systems  are  attracted  towards  the  same  dynamics 
even  though  the  initial  conditions  were  different.  The  predicted  dynamics  were  low-dimensional  in  nature, 
as  could  be  expected  from  the  severe  truncation  of  the  higher  modes,  and  tended  to  be  drawn  towards 
a  sinusoidal  limit  cycle.  The  dominant  frequency  appeared  around  500  Hz  which  is  low  compared  to 
measurements  performed  with  a  laser  doppler  anemometer  (LDA)  in  the  same  facility  at  x/D  =  6.  Hence, 
by  extrapolation  it  would  seem  natural  to  find  dominant  frequency  content  in  the  vicinity  of  1000  Hz  at 
x/D= 8,  the  present  location  under  investigation.  It  can  be  concluded  that  severely  truncated  cubic  models 
are  not  able  to  retrieve  the  higher  frequency  content  in  high-speed  flows  where  the  turbulence  spectrum  is 
broad  and  higher  modes  seem  to  play  an  important  role. 

The  next  four  LODS  are  developed  using  the  quadratic  form  of  the  ODE.  The  problem  being  much  smaller 
compared  to  the  cubic  model,  it  is  now  possible  to  include  many  more  modes  in  the  system  identification 
procedure.  This  is  done  with  the  aim  of  including  additional  smaller  scales  that  arc  truncated  in  the  cubic 
model  and  that  seem  to  be  crucial  in  highly  turbulent  flows  to  capture  the  higher  frequency  dynamics  which 
are  dominant  in  this  flow.  The  quadratic  model  is  able  to  be  contained  for  a  quite  long  period  of  time  (25 
ms)  in  a  bounded  domain.  Although  exhibiting  significant  periodicity  in  its  dynamics,  the  frequency  content 
is  shown  to  be  much  more  representative  of  the  actual  dynamics  existing  in  such  a  flow,  at  around  1000  Hz, 
intermittently  perturbed  by  higher  frequency  and  amplitude  fluctuations.  The  presence  of  high  frequency 
dynamics  arc  a  large  improvement  on  the  previous  cubic  systems,  which  is  largely  due  to  the  fact  that  a 
greater  number  of  modes  are  kept  in  the  modeling. 

LODS  4,  built  with  48  POD  modes  and  a  greater  number  of  training  samples,  is  not  able  to  stay  bounded 
for  as  long  as  LODS  3.  However  much  higher  frequency  dynamics  than  LODS  1  and  LODS  2  is  exhibited 
and  the  similar  range  of  frequencies  as  LODS  3,  around  1  kHz.  Most  likely  a  result  of  the  amount  of  TKE 
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truncated  from  the  system  in  the  identification  of  the  ODE  coefficients.  Purposely,  LODS  5  is  prescribed  a 
high  artificial  viscosity  of  to  study  the  consistency  of  the  effect  of  this  ad  hoc  parameter.  Indeed  the  POD 
coefficients  quickly  converge  to  constant  values  due  to  an  excessive  action  of  the  artificial  viscosity  in  the 
system.  The  process  of  finding  the  right  viscosity  to  compensate  for  the  truncation  is  quite  tedious  but  it 
is  indispensable  to  prevent  the  system  from  cither  diverging  (too  low  viscosity)  or  converging  towards  zero 
(too  high  viscosity).  However  once  the  parameters  of  the  dynamical  system  are  fixed  and  the  coefficients 
of  the  system  arc  identified,  the  successfully  tuned  viscosity  will  be  efficient  for  all  initial  conditions  and  at 
least  short  time  prediction  is  ensured.  LODS  6  is  developed  with  65  modes,  which  represents  the  maximum 
size  of  the  quadratic  system  without  it  being  ill-conditioned.  The  higher  number  of  training  samples  enables 
a  better  convergence  of  the  statistics  and  hence  a  better  conditioned  matrix  problem.  This  LODS  presents 
the  highest  frequencies,  which  is  consistent  with  the  finding  that  the  more  modes  arc  included  in  the  system, 
the  more  realistic  the  frequency  content  retrieved  in  the  system.  Additionally,  very  important  intcrmittency 
is  observed  in  this  case  with  high  amplitude  fluctuations  followed  by  periods  of  rather  low  fluctuations.  The 
study  of  the  frequency  content  of  the  predicted  velocity  field  demonstrates  that  a  severe  truncation  of  the 
higher  modes  in  the  case  of  the  cubic  modeling  prevents  any  high  frequency  dynamics  from  being  retrieved 
in  the  LODS. 

D.  Conclusions  and  recommendations 

The  moments  method  for  systems  identification  has  shown  that  it  can  be  applied  to  low  as  well  as  high 
Reynolds  number  flows.  The  only  shortcoming  of  low  order  modeling  at  high  Reynolds  numbers,  and  this  is 
not  specific  to  this  method,  is  that  the  smaller  scales  that  are  truncated  in  the  modeling  play  a  considerably 
more  important  role  for  the  dynamics  of  the  flow  as  in  the  case  of  low  Reynolds  number  flows.  Indeed  they 
arc  responsible  for  dissipating  the  TKE  produced  at  the  larger  scales.  Since  production  of  TKE  is  much 
larger  in  highly  turbulent  flows,  the  dissipative  scales  become  much  more  important  and  truncating  them 
presents  challenging  consequences.  This  is  not  as  much  an  issue  in  lower  Reynolds  number  flows  where  for 
one  part  a  lower  number  of  modes  captures  a  larger  part  of  the  TKE  and  also  where  the  dissipative  scales  arc 
not  as  crucial  since  production  is  not  as  important.  Hence,  in  the  low  Reynolds  number  flows,  a  cubic  model 
could  be  more  effective  at  predicting  the  dynamics  than  a  quadratic  one.  It  is  however  shown  here  that 
modeling  the  interaction  of  the  turbulence  with  the  mean  flow  (through  the  cubic  terms)  is  not  as  crucial  as 
using  as  many  modes  as  possible  (hence  140  as  much  TKE  as  possible)  in  the  development  of  the  dynamical 
systems.  It  is  therefore  recommended  that  quadratic  models  with  a  larger  number  of  modes  be  used  when 
developing  LODS’s  in  high-speed  and  highly  turbulent  flows. 

The  choice  in  the  truncation  level  has  shown  to  be  a  crucial  parameter  and  a  careful  study  of  the  spatial 
distribution  of  the  POD  modes  is  highly  recommended  to  prevent  the  system  from  being  unstable.  It  was 
shown  that  the  spatial  POD  modes  can  be  associated  to  others  by  similarity  of  spatial  features.  This  finding 
can  guide  the  choice  in  the  level  of  truncation  to  prevent  groups  of  modes  from  being  parted.  This  leads  to 
the  following  recommendation,  which  is  specific  to  axisymmetric  flows. 

When  the  measurement  permits,  it  is  suggested  not  to  assume  axisymmetry  a  priori  and  to  learn  from 
applying  the  2D-POD  to  all  three  components  of  velocity  instead  of  applying  the  Fourier  transform  in  the 
azimuthal  direction  followed  by  the  POD  in  the  radial  direction.  In  that  case  the  phase  relationship  between 
the  components  is  lost  and  the  interpretation  of  the  spatial  POD  modes  becomes  incomplete  and  perhaps 
misleading.  Indeed  interpreting  the  spatial  structure  of  the  POD  modes  has  shown  to  be  important  in  the 
choice  of  truncation  level  as  well  as  insightful  for  the  physical  interpretation  of  the  dominant  structures  in  the 
flow  in  terms  of  turbulent  kinetic  energy.  PIV  measurements  as  well  as  CFD  simulations  enable  estimations 
of  the  entire  field,  it  can  therefore  be  enriching  to  study  the  2D-POD  results  without  assuming  axisymmetry 
in  at  least  a  first  step. 
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